NO. SMI 
JANUARY 1956 


JOURNAL of the 


Soil Mechanics 
and Foundations 
VISION 


PROCEEDINGS Of THE 


AMERICAN SOCIETY 


OF CIVIL ENGINEERS 


VOLUME 82 


1 
‘ 
| 
il 
| 
— 
SORIETY OF 
cm 
ENGINEERS 
ql 


THIS JOURNAL 
... represents an effort by the Society to deliver technical data direct 
from the authors to the reader with the greatest possible speed. To 


this end, it has none of the usual editing required in more formal 
publications, 


Readers are invited to submit discussion applying to current papers. 
For papers in this iournal the final date on which a discussion should 
reach the Manager of Technical Publications appears in the table of 
contents and as a footnote with each paper. 


Those who are planning papers or discussions for “Proceedings” will 
expedite Division and Committee action measurably by first studying 
“Publication Procedure for Technical Paper” (Proceedings Papez 290). 
For free copies of this paper--describing style, content, and format— 
address the Manager, Technical Publications, ASCE. 


Reprints from this Journal may be made on condition that the full title 
of the paper, name of the author, page reference (or paper number), 
and date of publication by the Society are given. Copyright 1956 by 
the American Society of Civil Engineers. The Society is not responsible 
for any statement made or opinion expressed in its publications. 


This Journal is published quarterly by the American Society of Civil 
Engineers. Publication office at 2500 South State Street, Ann Arbor, 
Michigan. Editorial and General Offices at 33 West 39 Street, New 
York 18, New York. $4.00 of a member’s dues are applied as a subscrip- 
tion to this Journal. Application for second-class mail privileges is 
pending at Ann Arbor, Michigan. 


‘ 
3 
poe 
| 
| 
| 
| 
| al 
a 
oy ae 
f 
al 


JOURNAL 
SOIL MECHANICS AND FOUNDATIONS DIVISION 
Proceedings of the American Society of Civil Engineers 


CONTENTS 


January, 1956 
Papers 
(Discussion open until May 1, 1956) 
Paper 


861 Rupture Surfaces in Sand under Oblique Loads 
by Alfreds R. Jumikis eee @ 


862 Basic Concepts on the Compaction of Soil 


864 Control of Underseepage, Mississippi River Levees, 
St. Louis District, CE 
by C.I. Mansur and R. I. Kaufman .......... 


Structural Rigidity in Calculating Settlements 
by Samuel Chamecki 


Penetration Tests and Bearing Capacity of 


Cohesionless Soils, 


Triaxial Shear Tests on Pervious Gravelly Soils 
by Wesley G. Holtz and Harold J. Gibbs ... 


Division Activities 
1956-1 Newsletter 


1956-3 Newsletier ever eee eee 


| 
mm 4x. 
Page 
| 
861-1 
864-1 
| 865-1 
866 
867-1 
| 


Nag | 
q 
| | 
| 
¥ 
is 
“ae 
ae 
| 
Gis 
: 
ats 
1 
5 
| 
f 
‘Rie 
| 
| 


Paper 861 SM 1 


JOURNAL 
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RUPTURE SURFACES IN SAND UNDER OBLIQUE LOADS 


Alfreds R. Jumikis,! M. ASCE 
(Proc. Paper 861) 


SYNOPSIS 


The purpose of this article is to report and make available to the engineer- P. 
ing profession some of the results obtained by the author from experimental 3 - 
research in psommamechanics during the years 1939 to 1942. This research oy. 
pertains to deformation in soil, that is to say, the shape of the shearing or ens 
rupture (= sliding) surface in dry sand soil caused by an obliquely loaded a 
foundation model. This research shows that under the action of an excess Z 
oblique load there occurs a one-sided expulsion of sand mass from underneath f 
the foundation model, viz., the shearing off of a sand wedge. The sand wedge rs 
is definitely and clearly cut. It forms on the side of the eccentricity and in ie Pa 
the direction of action of the horizontal component of the obliquely applied na 
load. This is in contradistinction to a vertically eccentrically loaded founda- ie e 
tion model with no horizontal force component, where the sheared-off soil kc h 
wedge is formed on the opposite side of the eccentricity. Leg 
The shape of the observed shearing surface along which the one-sided ex- { 
pulsion or sliding out of the sand wedge takes place is a curved one. The in- 
tersection of the curved shearing or sliding surface with a vertical plane pro- 
jects a particular curve for the sand used in research. The observed 
particular curve of rupture agrees remarkably well with the mathematical b 
curve of a logarithmic spiral. The general equation of the spiral is given in i 7 | 
the text of this article. 
The sand wedge, sliding over the spiraled shear or sliding surface, re- 
sembles and acts as a solid body. a oo: 
This research presents certain features of interest to research workers in E: 


soil mechanics and to foundation engineers. It is also believed that it contrib- * 7 
utes to knowledge as to how and what kind of deformations take place in sand i 
soil under obliquely loaded foundations. It shows that the shearing surfaces 

are of a particular type of logarithmic spiral, depending upon the type of soil. * 7 


In addition, the knowledge of the shape of the sliding surface permits the engi- 
neer to perform stability analyses of obliquely loaded foundations, or to deter- 
mine the ultimate bearing capacity of soil under foundations loaded in this way. 


Note: Discussion open until May 1, 1956. Paper 861 is part of the copyrighted Journal ; — 
of the Soil Mechanics and Foundations Division of the American Society of Civil Engi- a 
neers, Vol. 82 No. SM 1, January, 1956. n 

1. Associate Prof. of Civ. Eng., Rutgers Univ., College of Eng., New Bruns- ee 
wick, N. J. 
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By knowing the shape of the sliding surface, the weight of the soil wedge can 
be taken into consideration in stability calculations. The properties of the 
logarithmic spiral likewise release the engineer from concern as to what law 
is followed by the reactions distributed underneath the soil wedge along the 
rupture plane. 

The observed shapes of rupture surfaces are valid for cohesionless or 
W-soil and for frictional-cohesive, y, c-soil as well. For pure cohesive or 
c-soil the spiral transforms mathematically into a circle. 


INTRODUCTION 


Methods of Loading . 

The evaluation of the bearing capacity of a soil along with its associated 
problem of stability and/or depth of a foundation is considered to be one of 
primary importance in soil mechanics and foundation engineering. 

In treating these problems of foundations supported on a level ground sur- 
face, the following methods of loading the foundation footing are usually en- 
countered: 

a) vertically and concentrically, 

b) vertically and eccentrically, and 

c) obliquely. 

It is known that one of the methods for the determination of the bearing 
capacity of soils and for stability calculations of the integral foundation sys- 
tem, “footing-soil,” is based on failure consideration, and by introducing a 
so-called “coefficient of safety” to avoid failure of the soil by shear. 

For this purpose, it is necessary to know, besides the physical properties 
of the soil, how the shear failure or rupture surface in the foundation- 
supporting soil takes place and what is the shape of the rupture or sliding 
surface. 

To learn this, research workers and engineers have suggested that in such 
calculations assumed shapes of failure surfaces be used. It is the practice to 
use these assumed shapes rather than to determine them analytically, or to 
base calculations on experimentally observed rupture surfaces. The reason 
for this is that it is difficult to determine the direction of the principal 
stresses in a fragmental medium such as a mass of soil. Because soil is an 
indeterminate and therefore a difficult material to work with, research work- 
ers usually performed their investigations relative to this subject on a small 
scale and with a cohesionless material, namely, dry sand. This permits ex- 
cluding from their studies the effects of moisture and apparent or transitional 
cohesion associated with it. A further factor to be considered is moisture 
and its migration with temperature changes. 


The Shape of Rupture Surfaces in Sand under Various Methods of Loading 


Rupture Surfaces Under Vertically and Concentrically Loaded Footings 

The problem of the shape of the rupture surface in dry sand under a verti- 
cally concentrically loaded footing or loading plate is frequently dealt with in 
the special technical literature. However, it must be said that even now this 
problem has not yet been completely solved. The rupture or sliding surface 
in the soil resulting from failure of soil by shear was thought of by various 
authors as having several differing shapes, each according to his own theory. 
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Only during the past half century has the problem of stability of vertically and 
concentrically loaded footings been studied more seriously. The question is 
theoretically discussed by Résal and Frontard, (13) Schwedler and Zimmer- 
mann,(39) Prandtl,(28) Pihera, (27) Fréhlich, (10,11) Rendulic(29) and others, 
. whose formulas are also mentioned by Terzaghi. (34,35) Prandtl’s theory as 
applied te soil mechanics is described also by Krynine,(16) Taylor (33) and 
Tschebotarioff.(38) The problem of small scale bearing capacity studies un- 
der vertically and concentrically applied footing loads has been treated ex- o yt 
perimentally by Kurdumoff,(17) Krey,(15) Fellenius,(7,8) Housel, (12) a 
Meischeider,(19) Jumikis,(4) M, G. Bekker, (1) Peynircio$lu, (25) Caquot and 
Kérisel,(2) Davis and Woodward, (5) and possibly others. 
Specific shapes of the rupture surface were assumed by the following 
authors to be these: 
straight (plane) by Pauker,(6) Rankine, (34,35) Ritter; (32) 
i circular by Fellenius,(7,8) Krey,(15) Chatley; (3) = 
circular or combined with a tangent by Krey, Terzaghi; (15) oF 
cycloidal by Résal,(30) Frontard, (9) Collin; (4) 
logarithmic spiral by Schwedler, Zimmermann, (39) Frohlich, (10) 
Pihera,(27) Jumikis,(14) McLeod, (18) Davis; (5) ue 
combined, that is, logarithmic spiral with a tangent, by Prandtl, (28) Caquot oo 
and Kérisel, (25 Résal, (30) Reissner. (31) 


on Rupture Surfaces Under Vertically and Eccentrically Loaded Footings 
The shape of the rupture surface in dry sand caused by vertically and Re 

eccentrically applied loads was studied experimentally by Jumikis in 1939. (14) a 
+ Its shape was found to be a logarithmic cylindrical surface. 4 rs 
The bearing capacity of foundations under eccentric loads was studied ex- i 
. perimentally and theoretically (in 1948 and 1953) by Meyerhof.(20,21,22) The = 
curve of rupture surface used in his studies is spiraled. A note by Meyerhof 7 
(22, page 441) indicates that Ramelot and Vandeperre did experimental re- 

search in Belgium in 1950 on eccentrically loaded foundations. 


Rupture Surfaces Under Obliquely Loaded Footings 

According to the author’s best belief, there is really very little technical = 
ia information available or accessible relative to the process of deformation = } = 
hed which takes place in soil under obliquely loaded foundations, as well as the |; 
a shape of the rupture surface. Before 1944, as far as the author knows, there 
. are probably only three authors whose works have some bearing upon this 
question. They are Moller, (24) who described the measurements of the re- 
me sistance of sliding of construction materials on the soil in connection with 
making decisions for the Cologne-Muhlheim bridge system across the Rhine 
River; Pialoux,(26) who measured the shearing resistance of the soil in con- <i 
nection with building railroads in the French colonies, and Terzaghi,(37,38) ft 
who determined the resistance to sliding of an anchor block of a suspension ae : 
bridge. 

However, every civil engineer encounters obliquely loaded foundations in 
his daily work, for example, wherever in addition to vertical load, norizontal Bt 
forces also act on a foundation. Such cases are foundations of vaults, arch- 

. type bridges without ties, anchor blocks of suspension bridges, retaining walls, 
of frameworks with horizontal reaction-components, dams and hydraulic struc- 

tures. 

In order to perform calculations for the bearing capacity of a soil, for the : 

7 depth of a foundation or over-all stability of a foundation acted upon by oblique F 
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loads, it is necessary to know the mode of failure in shear in the soil under- 
neath the foundation. 
Because the problem of the stability of obliquely loaded foundations was 
and still is insufficiently studied, in 1939 the author undertook in the labora- 
tory experimental research with small scale models supported on dry sand 
and loaded obliquely.(14) The observed rupture surfaces were found to be 
logarithmic spirals. 
The small scale experimental apprcach of study was chosen because 
1) it is easier to handle the testing apparatus; 
2) it is less expensive to procure testing equipment and to perform the ex- 
periments, than if this work had to be executed to a large or natural 
scale; 
3) it is possible to eliminate many other factors involved, and to control 
and observe only the ones in which we are particularly interested; 
4) by artificially simulating field conditions in the laboratory 
a) it is possible conveniently to observe the processes and shapes of 
defor mations in prototype soil, and 

b) it is possible to perform actual measurements of deformations in the 
soil prototype under given loading conditions, thus obtaining the data 
for the correlation of pertinent facts; 

5) the researcher, in performing stability calculations of foundations, is 
released from making several assumptions usually made in the mathe- 
matical approach of research (for instance, that soil is a perfect elastic 
material with isotropic properties, which it is not). 

At this point must be mentioned the experimental-theoretical works on 
preety, capacity of foundations under eccentric and inclined loads by Meyerhof 
(1953), 22) his work on foundation research and its application to design 
ose as well as his work on bearing capacity of shallow footings on dry 
sand. 


Object of Research © 


Definition 
The definition of the author’s research was formulated as follows: 
“Experimental research concerning the shape of the rupture surface 
in dry sand caused by an obliquely loaded foundation model.” 


Objective 
Previous experimental research by the author with vertically and eccen- 
trically loaded models showed that when the shear strength of the sand used 
for experimentation was exceeded, a sand wedge underneath the footing was 
sheared off along a curved failure surface. This was expelled from under- 
neath the footing as a solid body on the opposite side of the eccentricity. It 
was, therefore, the object of this research to find out 
1) whether under the action of a horizontal force component on the footing 
model a soil wedge can be expected to be pushed out from underneath 
the footing, or whether, because of the expected lateral sliding of the 
footing model over the ground surface, there is no soil wedge in the 
soil sheared off at all; 
2) how, if a soil wedge were sheared off in the soil mass, supporting the 
footing model, this process takes place; 
3) what would be the shape of the rupture, viz., sliding surface, and 
4) what is the dependence of the sliding surface upon the factors which 
partake in the experiments. 
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Paper 861 JUMIKIS 
Methodology 


Apparatus 
The apparatus for carrying out the experimental research consists of 
1) a 30" x 6.5"'x 12" large steel box with glass walls for filling and loading 
the sand, 
2) a concrete foundation model, 4" x 6'"'x 6" in size, and 
3) a loading device for applying to the model vertical and horizontal loads 
by means of dead weights (Fig. 1). 


Sand 

For tre reasons mentioned in the introduction, the part of the experimental 
research here described was carried out with dry sand. It was a typical 
Danube sand, the grain size distribution of which is represented in Fig. 2. 
The specific gravity of this sand was Y, = 2.67, the coefficient of its internal 
friction was found to be tanw = 0.580, and the coefficient of friction between 
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the concrete and the sand at the base of the concrete foundation model was 
tany, = 0.630. 


Method of Detecting Rupture Surfaces 
For the detection of the rupture surface in dry sand caused by an obliquely 
loaded foundation model, black-stained horizontal sand layers, 3 mm thick and 
spaced vertically 10 mm apart, were tamped in the box. An obiique load was 
effected by a constant vertical load and a progressively increased horizontal 
load (Fig. 3) until failure by shear in sand occurred. At this point a one-sided 
sheared-off sand wedge underneath the base of the model and “ground surface” 
showed up nicely in the form of fault- or shear-lines. These appeared in the 
form of a displacement of the black, horizontally placed lines, within the 
wedge as compared with those in the rest of the sand mass, indicating that the 
sand wedge rotates along the cylindrical rupture surfaces. 


Vv 


FOUNDATION MODEL, 4°X6"x6" 
SAND BOX, 


PRINCIPLE OF OBLIQUE LOADING. 


FIG. 3 


These rupture surfaces were photographed, analyzed, mathematically and 
grapho-analytically treated, and polar equations for these rupture surfaces 
were determined. 

After the photographs and other data were taken, the model of the founda- 
tion and the sand were removed from the box. The box was then refilled with 
sand at the proper density, the “ground surface” smoothed, the model put into. 
position and again subjected to loads. 

The applied vertical loads exerted 0.25, 0.50, 0.75, 1.00 and 1.25 kg/cm? 
(t/ft?) contact pressures on the sand. At each of these contact pressures the 
model was subjected to horizontal loads at the following height above the base 
of the model: 

h = 0, 2", 4" and 6"", respectively. 

The method of loading the foundation model in most cases corresponds ap- 
proximately to that encountered actually in engineering practice, namely, first 
the weight of the structure increases as the structure goes up, then the 
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horizontal loads (lateral earth pressure, hydrostatic pressures, wind loads, 
braking action of rolling stock on railroads, ship impacts on quays and others) 
gradually begin to apply. 

Upon the gradual application of a horizontal load to a vertically and con- 
centrically loaded foundation model the latter displaces in the direction of the 
horizontal load, and continues to settle. Upon further increment of the hori- 
zontal load the model tilts by a clockwise (in this case) rotation, creating a 
conventional trapezoidal pressure distribution with the maximum pressure 
under the edge of the model on the side of the eccentricity. The sand on this 
latter side is now denser than before, and upon exceeding its shear strength 
a soil wedge, resembling a solid body, is sheared off, rotates and slides out 
along the shearing surface from underneath the model. The shearing of the 
sand wedge takes place suddenly and is accompanied by a slithering noise 
characteristic of sand. Figs. 4, 5 and 6 illustrate the principal difference be- 
tween the shearing deformation in dry sand of a vertically-eccentrically 
loaded foundation and that which appears under an obliquely-eccentrically 
loaded foundation. 

Figs. 7, 8 and 9 show rupture surfaces under oblique load at contact pres- 
sures of 0.50, 1.00 and 1.25 t/ft2, respectively. The point of application, h, 
of the horizontal load in this series of experiments is zero. 

From these figures it can be seen where the sliding surface starts, how 
deep it reaches below the ground surface, the nature of its curvature, and 
what is its extent laterally. Also, as can be observed, the outermost branch 
of the sliding curve near the ground surface, beginning from a certain point, 
is, in some instances, bent somewhat down in a position M—N, making a 
“creased” impression (Fig. 10). This is partly explained by the unconfined 
boundary condition at the ground surface. 

Thus, the method of detecting the shape of the rupture surface in sand is 
based on static loading of the foundation model and on actual observations of 
the behavior of the foundation model supported directly by a soil, nor does it 
require any theoretical assumptions as to its probable behavior. 


The Shape of the Rupture Surface 


Analysis of Spirals 

Examination of the observed shapes of rupture surfaces indicated that the 
curve which forms intersecting the sliding surface with a vertical plane could 
be most easily expressed in a polar coordinate system. This proved to be 
true comparing the observed curves with theoretically calculated ones, using 
some common radius-vectors, r, and amplitudes,w. It revealed that the rup- 
ture surface curve coincides remarkably well with a curve of a logarithmic 
spiral, the general equation of which is 


(la) 


(1b) 


where 
r = radius~vector, 
lg = reference vector, or a segment on the polar axis from the pole of 
the spiral cut off by the spiral atw= 0, 
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e = base of natural logarithms, 
w = amplitude, or angle between r,, and r, and 


tany = tan( 5 -v) 
tanw = coefficient of internal friction of sand. For the sand used in 
this research tan = 0.580; 


(z -y )- angle between a radius-vector and its corresponding normal 


of the spiral. 
The minus sign at w denotes that as the amplitude w increases, the radius- 
vector decreases (Fig. 11). This method was adapted for the convenience of 
analyses of experimental data. 
Analysis of two mathematical logarithmic spirals with the common tan ~ 
indicates that 


Trane const, (2) 


and ¥. (4) 


where ry is the radius-vector when w= 7. 


O - Pole. OX - Polor axis. m = tan ton . 
Radius - Vector. w Amplitude . 


SPIRAL . 
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In such a case (eq. 2, 3, 4) two spirals are similar and therefore the center 
of similarity is the pole of the spiral (Fig. 12). 

These and other properties of the logarithmic spiral were used for studying 
the observed rupture curves, for determining the position of poles of the ex- 
perimentally found spirals, and for setting up the physical equation for the ob- 
served spiraled rupture curve. 


= const. 


FIG. 12 SIMILAR SPIRALS 


Comparison of Theoretical and Observed Spirals 

After determining the positions (coordinates) of the spiral poles, the 
lengths of the radii-vectors at -w = 15° increments were computed theoretical- 
ly, and measured on experimentally obtained spirals (Figs. 13, 14). These 
radii-vectors were compared and tabulated. For example, Table 1 shows 
radii-vectors which were found experimentally (rg) and theoretically calculat- 
ed (rc) for five spirals resulting from contact pressures of 0.25, 0.50, 0.75, 
1.00 and 1.25 t/ft2 for amplitudes of 15° increment. In this instance the hori- 
zontal load, H, was applied h = 4" above the base of the foundation model. 

A procedure similar to that just described was also applied to models 
acted upon by a horizontal load at different heights of application other than 
h= 4" 

A comparison of the lengths of the theoretically calculated and experi- 
mentally obtained radii-vectors is shown in columns 4, 7, 10, 13 and 16 of 
Table 1 by way of differences and in percent. As can be seen from the photo- 
graphs, the values of r150, T1g5 and rygo are difficult to measure. The com- 
parison ot calculated and observed radii-vectors (r, and rg, respectively) 
values, except for the three for - w = 150°, 165° and 180°, are in close agree- 
ment, which permits the conclusion that the observed rupture surfaces coin- 
cide practically very well with logarithmic spiral curves. It can be observed 
that the agreement of the r, and r, values depends very much upon the ac- 
curacy of locating the position of the pole of the spiral. 
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Similarity of Spirals 

Based on the r¢ values (Table 1), the similarity of the five spirals in ques- 
tion is shown in Fig. 15. Evaluation of data shows that for the experimental 
series with the applied five contact pressures and for the applied horizontal 
loads at h = 0", 2"'and 4" above the base of the model, and for the purpose of 
the analyses, the important radii-vector ratios are as follows: 


Here rgg indicates that radius-vector the amplitude of which is -w = 60°. 

Hence, for analyses and for establishment of the physical equation of the 
rupture surface for the sand used, the following average ratio-values can be 
used: = 1.82, and = 6.20. 


The Size of the Logarithmic Spiral 

From equation (1) it can be seen that the size of the spiral depends upon 
tany, (-), and Tq: But experiments show that the size of the spiral depends 
also upon vertical and horizontal loads, V and H, viz., their resultant R, act- 
ing upon the foundation model, and the width of the latter, as well as upon the 
eccentricity of V and the inclination of R. The effect of the width of the found- 
ation model will not be discussed here. The influence of the eccentricity and 
the inclination is included in the loading results. However, they are masked 
out in this presentation. Because the measurement of r, from the rupture 
surfaces is uncertain, a radius-vector which could be measured with a greater 
accuracy than ro is used for the determination of r,. Since the rupture sur- 
face shows up most clearly in the central part of the curve, the radius-vector 
rgg (at -w = 60°) was chosen for determining ro, as well as for other studies 
of the observed spirals. 

Analysis revealed that for the experimental conditions the ratio of the 
reference radius-vector to that at an amplitude of -w = 60° is constant, 
namely: r,/T60 = 1.82 = constant, which gives the desired rp: 


Ty = (1.82) - (5) 


Further, it was found that the general equation of rgo in terms of the resultant 
load, R, is: 


19 = 
in r r 
: 60 TT 
0 1.82 6.20 
211.82} 6.20 
4 | 1.82] 6.28 
| 
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Substitution of eq. (6) iato eq. (5) gives: 


a 


ry = (1.82) k +R (7) 


The quantities k and a are constants and depend upon the sand, width of foot- 
ing and the loading conditions. The quantity R is the resultant of V and H: 


Weer. (8) 
Equation of Spirals 


The general equation of the various logarithmic spirals observed in this 
research can be expressed by the following formula: 


r = (1,82) > k- R& tan¥, (9) 


For h = 4", k = 6.31 and a = 0.47, and with tany = 0.58, the particular equation 
of the logarithmic spiral is 


Hence, this research shows that the shape of the rupture surface in dry sand 
caused by an obliquely loaded foundation model can be considered as being a 
particular and uninterrupted (continuous) logarithmic spiral for the sand used. 

Also, the shape of the rupture surfaces in sand under a vertically- 
eccentrically loaded foundation model is a particular logarithmic spiral for 
the sand used, which differs somewhat from that as expressed by eq. (10). 


CONCLUSIONS 


The described experimental research shows that 

1) from underneath an obliquely loaded foundation model a sand wedge is 
sheared off and expelled as a solid body laterally one-sidedly along a 
rupture surface; 

2) the shearing off of the sand wedge upon exhaustion of its shearing 
strength takes place suddenly; 

3) the shape of the rupture, viz., sliding surface in dry sand, caused by 
vertically-eccentrically and obliquely loaded foundation models is a 
cylindrical surface, of a logarithmic spiral shape; 

4) the general polar equation of the logarithmic spiral curve for an 
obliquely loaded foundation model is 


r = (1.82) - k + R® - @-W> tang, 


: 
; 
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5) the knowledge of the shape of the rupture surface permits the engineer 
to determine the bearing capacity of the sand soil, or to calculate the 
depth (viz., surcharge) for laying of foundations in sand soil, or to per- 
form stability analyses of footings on sand soil; 

6) with various resultants, R, and with different heights, h, of application 
of horizontal loads, H, the spirals, representing rupture surfaces are 
similar; 

7) under obliquely loaded foundation models the poles of the spirals are 
located above the “ground surface,” and in the front third (direction of 
motion) of the width of the base of the footing; 

8) at failure, under oblique loading, the entire footing or foundation model 
does not rest on its entire width on the sheared-off spiraled sand wedge 
(Figs. 7, 8, 9); 

9) by a vertically-eccentrically loaded foundation model the starting point 
of the rupture surface in sand is practically situated at the edge of the 
model on the side of the eccentricity, and the model rests completely 
on the sliding sand wedge; 

10) the sizes of the spirals depend upon the R- magnitudes, namely, rgg- 
radii are directly proportional to the resultant loads (viz., functions o 
V and H); 

11) the horizontal load H, necessary to produce a rupture surface, depends 
upon the vertical load, V, as well as upon the height of application, h, 
of the horizontal load; 

12) the lower the height of application, h, of the horizontal load, H, is, and 
the larger V is, the more clearly pronounced are the rupture surfaces; 
at a certain high value of h, the overturning phenomenon of the medel 
starts and no rupture surface can be observed. Here, for each vertical 
load, V, it is possible to find height, h, at which a horizontal load of 
very small magnitude (nearly zero) is sufficient to cause a condition 
impending to tilt; 

13) in the case where the sand soil contains also some cohesion, it is pos- 
sible by means of a simple integral to perform stability calculations 
(see Fig. 16). The moment of cohesion relative to the pole of the spiral 
is: 


M =—%—(r 2) ; (11) 


14) theoretical considerations indicate that when tany = 0, i.e., if we have 
to deal with a pure cohesive soil, this spiral equation is transformed 
into a circle: 


r = = const; (12) 


15) by means of similarity or model law, the logarithmically spiraled rup- 
ture surface can also be used for stability calculations of large scale 
foundations. 

16) There are no difficulties in drawing spirals, particularly when spiral 
tables or spiral templets are applied, or the rules for the construction 
of spirals are employed. 
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of internal Friction. Ry Reaction. 
Gohesion: c, in pounds per unit srea; 


Cohesion : (C)(ds)(i) (pounds per differential area). 
Moment : dM,* c-ds-cosP-r=c- dw. 


FIG. 16 MOMENT OF COHESION 


The application of a logarithmically spiraled rupture surface to stability 
calculations possesses a great advantage over other assumed surfaces, name- 
ly, that it is not necessary to assume nor to be concerned as to how the reac- 
tions in sand soil along the sliding surface are distributed and what their 
magnitudes are. This is because the logarithmic spiral possesses an im- 
portant property, namely—all radii-vectors pass through the pole of the spiral 
(moment arm for reactions is 0). Hence all soil reaction moments are auto- 
matically excluded from our stability calculations vy comparing active (driv- 
ing) and reactive (resisting) moments. 
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SYNOPSIS 


An attempt has been made to discuss the soil compaction process ona 
rational basis and to create basic concepts of the problems involved. The 
Law of Conservation of Energy is employed to interpret how energy is being 
spent during the process. The behavior of a soil mass during compaction due 
to the influence of various factors is explained. The “Stage Compaction” con- 
cept is introduced. The effect of gravel content on the compaction of fine soil 
fractions is presented. The practical sign‘ficance of this rational approach is 
also included. 


INTRODUCTION 


The available information on soil compaction pertains mainly to the per- 
formance of methods rather than to the fundamental principle of compaction. 
Little comprehensive soil compaction investigation has been made with a view 
to establishing a basic principle or theory governing the relationship among 
the numerous variables that influence the compaction result. Engineers facing 
the problem in construction must rely principally on past experience as basis 
to make decisions and therefore are subject to the shortcomings of all empiri- 
cal methods. Controversial and divergent views are frequently raised and ag- 
gravated on important subjects such as: the choice of types of roller, the 
proper use of them, the maximum soil density that can be attained economical- 
ly, the method of field control, the laboratory Proctor Test, the compaction of 
gravelly soil, etc. The significance of the problem cannot be overemphasized, 
especially the economical aspect, when the extent and scale of earthwork con- 
struction is realized. A fundamental understanding of the compaction problem 
and the formation of basic concepts are necessary in order to correctly pre- 
dict and efficiently control the influence of the various factors on the final re- 
sult. 


Note: Discussion open until May 1, 1956. Paper 862 is part of the copyrighted Journal 
of the Soil Mechanics and Foundations Division of the American Society of Civil Engi- 
neers, Vol. 82 No. SM 1, January, 1956. 

1. Civil Engineer, Gannett Fleming Corddry & Carpenter, Inc. Harrisburg, 
Pennsylvania. 
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Assumption 


The following study on compaction will be based on the assumption that the 
soil mass, in its loose state before compaction, is structurally homogeneous. 
It is free from large chunks, lumps or large bridged voids. It is assumed 
that before rolling, the soil is subjected to the usual procedures of excavat- 
ing, transporting, dumping and spreading in layers. The passing of hauling 
and spreading equipment on the loose soil may result in a reasonably struc- 
turally homogeneous soil. However, it must be recognized that a natural soil 
mass is never truly homogeneous either in composition or in structure. It is 
believed that the simplified assumption of homogeneity will introduce no spe- 
cial inconsistency when dealing with the general trend of the property or sta- 
tistical average results of the actual fill. 


Compaction Defined 


When a soil mass is subjected to a loading, whether static or dynamic, 
transfer stresses and subsequent strains are involved. The strains produce 
deformation and volume change. A definite period of time is required for the 
non-elastic strains to develop completely in a soil mass after the inducement 
of the stresses, due to the requirements of drainage and plastic readjustment. 
If the loading is applied long enough to overcome this time lag, the soil is 
said to be “consolidated” with respect to the load. However, if the loading 

is instantaneous, such as the blow of a hammer or passing of a roller, there 
is only time enough for the occurrence of initial strains or partial compres- 
sion. The latter process is called “compaction.” Generally, compaction may 
be defined as the process of densification of soil by dynamic load application, 
causing a decrease in air voids due to cuange in relative positions of soil 
grains. 


Law of Conservation of Energy 


The stress-strain velationship of a soil mass during compaction is an ex- 
ceedingly complex problem. It is improbable that it will ever be expressed 
by formulas which take into consideration all the major influencing factors. 
Nevertheless, the first law of thermodynamics, well known as the Law of 
Conservation of Energy, may be employed to explain qualitatively the basic 
factors involved during the compaction  prewves 

The Law of Conservation of Energy(1) states that “Energy can neither be 
created nor destroyed but only converted from one form to another; or, for 
any system, the net result of heat or work will be a change in the energy of 
the system.”2 Thus 


E,-€E, 


=AE 


(1) 


2. (The System may be defined as the region where transfer of mass and 
energy are to be studied.) 
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where Q = heat added to the system; 
W = work done by the system; 
E, and E, = initial energy and final energy of the system 
respectively; 
AE = change in total energy of the system. 


Equation (1) may be written in the form of 
(-W) = AE + (-Q) (2) 


which means that the work done on the system is equal to the change in total 
energy of the system plus the heat transferred from the system. 

In the process of com»action of soil, the terms of the above stated equation 
may be explained as follows: 

a) Work Done on the System - Work is defined as the product of force and 
the distance through which the force acts. In static loading, the work done on 
the soil is 


w={Pds (3) 


where P = load on the system; 
S = settlement. 
In dynamic loading, the work done is 


w=tM(ve- vi?) 


where M = mass of the applied load; 
v1 = initial velocity of the applied load; 
vg = final velocity of the applied load. 

Generally, the total amount of work done on the system can be determined 
or estimated. In laboratory soil compaction by the Proctor Method, the pro- 
duct of the hammer weight and the distance through which the hammer drops 
is the amount of work done, and is the energy transferred to the soil sample. 
In field compaction, the work done on the soil varies directly as the draw-bar 
pull on the roller and the distance through which the roller travels. 

b) Change in Total Energy of the System(2,3) - The energy of a system re- 
fers to the internal energy that the atoms and molecules of the system possess 
as a result of their configuration and motion. 

In soils containing particles in clay sizes, the increase in density by com- 
paction changes the total internal energy of the soil. It is a process involving 
reorientation of soil particles which possess forces of repulsion or attraction 
due to their adsorbed ions and adsorbed water molecules. The physical 
changes during compaction, such as the compression of gases in the voids, 
the increase of the amount of dissolved gases in pore water, the elastic strain 
of solid particles, and other complex colloidal phenomena, are the result of 
change in the amount of internal energy. Increased energy is being stored in 
the system but the process is reversible. The amount of the stored energy 
in cohesive soil is demonstrated by the rebound or swelling upon the removal 
of load. 
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The compaction of cohesionless soil is a much simpler physical process, 
in that the individual soil grains are forced closer together. The change in 
total energy of the system is principally the inducement of elastic strain 
energy. 

c) Heat Transferred from the System - In the process of compacting soil, 
there must be relative motion of soil particles of all sizes. Force is required 
to overcome the frictional resistance developed between particles during the 
motion. As is well known, the energy spent to overcome the frictional resist- 
ance is converted into heat, which dissipates into the surrounding media. 

There are three types of frictional resistance which characterize the rela- 
tive motion between soil particles during compaction, namely; (1) Dry surface 
friction; (2) Hydrodynamic or thick-film lubricated surface friction; and 
(3) Boundary or thin-film lubricated surface friction. 

1) Dry Surface Friction(4) - The frictional resistance between two dry sur- 
faces is determined by the coefficient of friction and the normal pressure. It 
is considered independent of area of contact and of relative velocity of motion. 
The compaction of dry granular soil consists of overcoming the dry surface 
frictional resistance between soil particles which arises from sliding, rolling 
or interlocking. 

2) Hydrodynamic or thick-film lubricated suriace friction(4) - When a 
lubricant such as oil is being gradually added to the dry surfaces of sliding, 

a lubricative film will be formed over the dry surfaces. It is very thin at 

first and increases in thickness as more and more lubricant is present. 
Hydrodynamic lubricated surface friction will occur when there is a plentiful 
supply of lubricant to form a film, the thickness of which is large in compari- 
son with the height of irregularities of the surfaces. The frictional resistance, 
which is independent of the materials composing the sliding surfaces and of 
the pressure between them, is entirely different from that of the dry friction. 
The few factors which determine the force required to produce the relative 
motion may be expressed by the following equation: 


Av 
(5) 


Where F = force required to cause the motion; 
HM = Coefficient of viscosity of lubricant; 
A = area of the surface of motion; 
Vv = relative velocity of the surface of motion; 
h = distance between the surfaces of motion. 

This type of surface condition exists only in very fine-grained soils such 
as Clay or fine silt because water, the common lubricant for soil, will not 
form a thick enough film on coarse grain. For compaction of clay soil, Equa- 
tion (5) may be viewed as representing the statistical average value of the 
hydrodynamic frictional resistance. It can thus be seen that the force re- 
quired to overcome the frictional resistance varies directly as: (1) the 
viscosity of water which is constant at a given temperature; (2) the surface 
areas of the particles; (The fact that one percent of clay in a sand contributes 
more than ninety percent of the total internal surfaces may indicate the in- 
fluence of clay content on the energy of compaction); and (3) the velocity of 
motion, which shows that more energy is required for dynamic compaction 
than that for static compaction. The force to cause motion is inversely pro- 
portional to the thickness of water layer between particles. 
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3) Boundary or thin-film lubricated surface friction - Between the extreme 
conditions for dry friction and hydrodynamic lubricated surface friction there 
is boundary or thin-film lubricated surface friction. It occurs when the supply 
of lubricant is not enough to reach the hydrodynamic surface condition. How- 
ever, the lubricant tends to reduce the coefficient of friction between the dry 
surfaces. The contact pressure and the nature of the materials composing 
the bodies are still of importance but the required force is relatively inde- 
pendent of the area of contact of surfaces and the velocity of slicing. This 
may be considered as the transition of the above discussed two extreme con- 
ditions and may apply to soil with intermediate grain size such as fine sand 
and silt. 


Significance of Law of Conservation of Energy 


The significance of the Law of Conservation of Energy upon compaction of 
soil lies in understanding how the energy is being spent during the compactive 
process and how it may be consumed in many other ways than the densifica- 
tion of soil. Two factors of primary importance are frictional resistance of 
soil, and unaccounted for or “wasted” energy. 

Since compaction is the result of reducing voids in a soil mass, it must 
involve the rearrangement of soil grains by shifting their relative positions. 
Energy is required for overcoming the frictional resistance of movement. 
However, the rearrangement of soil particles does not necessarily follow the 
reduction of air voids, due to the complete shearing failure of the soil mass 
to resist the load. The latter is one form of “wasted” energy. 

The repeated rolling of a light roller upon a dense soil may cause very 
little or insignificant decrease in air voids. However, in accordance with the 
Law of Conservation of Energy, none of the energy is lost; it occurs in the 
soil mass, similar to an elastic structure, as potential energy, usually desig- 
nated as strain energy. From the viewpoint of compaction, strain energy is 
entirely “wasted” energy. 

Other forms of “wasted” energy which are inherent during compaction are 
those for elastic strain of soil grains, for compressing air in voids, for over- 
coming bond among particles and for other possible physio-chemical reac- 
tions. 

It seems clear that only by achieving the reduction of both frictional resist- 
ance of compaction and “wasted” energy can the efficiency of the compactive 
process be improved. The practical approach to the problem will, of course, 
be the understanding of the important factors which influence the result of 
compaction. The following discussion of “Loading and Resistance” may form 
a basic concept to explain the influence of various factors upon compaction 
which will appear in subsequent sections. 


Loading and Resistances 


The modern large scale field compaction of soil is done generally by the 
rapid penetration and withdrawal of the tamping feet uf sheepsfoot rollers or 
of the highly inflated rubber tires of pneumatic tire rollers. This rapid ap- 
plication and removal of loading is essentially dynamic.3 Dynamic loading is 
quite different from static loading due to the fact that it produces so-called 


3. (Compaction by vibration will not be discussed in this paper.) 
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dynamic resistance in the soil in addition to the static resistance. 

The variation of resistance of soil during dynamic loading is shown in 
Fig. 1-A.(5) The dynamic resistance reaches a maximum under the impact 
of the applied load. It decreases rapidly and approaches static resistance 
and final equilibrium with the load at the end of settlement. The stress and 
strain distribution under the load may be regarded as similar to that for 
static loading. The pressure bulb concept, the effects of size and shape of 
loaded area upon the soil resistance and the plastic equilibrium theory may 
be employed to explair the behavior of soil mass during compaction. It should 
also be emphasized that the strains or settlements are a direct function of 
loading intensity rather than the gross load. 

If the applied stresses exceed the shearing strength or resistance of the 
soil, local failure begins and the load starts to sink into the soil. As the load 
sinks, the soil under the load is pushed downward and outward. This motion 
will mobilize more and more resistance, consisting not only of the increased 
resistance due to lateral confinement from depth but also that due to the in- 
crease in soil density that results from the settlement motion itself, provided 
the soil mass is not completely saturated. The settlement stops when equilib- 
rium between stresses and resistance is reached. 

If the resistance of the soil is relatively high compared with the stresses, 
the load will cause very little settlement, which also means very small amount 
of void change and, thus, little compaction. Most of the energy is spent in 
producing shearing and volumetric strains. Only a small amount is used to 
attain compaction through the movement of particles. If, however, the re- 
sistance is extremely low in comparison with the stresses, the load will cause 
a complete shear failure of the underlying soil by sinking deep and fast and 
replacing the volume of soil by pushing it in an outward direction. This com- 
pletely disturbed state may result in compaction of the soil under the load and 
loosen the soil on the sides. The total net reduction in voids is questionable 
and is inefficient at best. Energy is spent in compacting one portion of soil 
and, at the same time, loosening another portion. Both compacting and loosen- 
ing involve movement of particles and require energy to overcome the fric- 
tional resistance. 

From the above discussions it may be seen that the compaction depends 
upon the relative values of load and resistance, and upon the amount of settle- 
ment. 

It can be said then that the optimum roller pressure seems to be that 
which is large enough to cause local failure without proceeding into complete 
failure of the soil mass. This has been borne out by actual experience. For 
example, the passage of the heavy sheepsfoot roller(6) in one case caused 
complete shear failure of the soil under the roller but produced very little 
compaction of the soil mass, but by reducing the roller pressure, satisfactory 
densities were obtained. 


Introducing “Stage Compaction” 


When the loose, cohesive or moderately cohesive soil is first spread out 
in layers ready for rolling, its shearing resistance is always very low. It is 
likely to be incapable of supporting a heavy roller without complete failure. 

In using the sheepsfoot roller, the tamping feet completely penetrate the un- 
compacted soil, and the roller drum comes in contact with the soil, effectively 
spreading the load to a large area. In using the pneumatic tire roller, the tire 
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sinks into the soil and reduces the loading pressure by increasing the contact 
area. It is seen that the soil will be compacted to some extent from the first 
passes cf the roller. Subsequent trips of the roller are thus on a soil with 
highvr resistance and every additional pass of the roller is on a stronger soil 
thaa that of the previous one with the result that it becomes possible to 
achieve a certain maximum density for each particular loading. 

It is thus seen that any of the familiar heavy rollers is likely to be too 
heavy at the beginning of rolling loose soil layers. However, the efficiency 
of the roller gradually increases until a maximum is reached. Then, the ef- 
ficiency decreases because the roller pressure is too light relative to the re- 
sistance of the soil. At last, the efficiency of roliing approaches zero because 
no further compaction will result from additional rolling. That is, the shear- 
ing strength of the soil is equal to the roller pressure. The work done by 
- rolling is stored entirely as elastic strains. Fig. 2-A shows the rate of 
density increase with respect to the number of roller passes. It is a measure- 
ment of the efficiency of rolling. 

Light pressure rollers have a much higher starting efficiency than the 
heavy pressure rollers. However, a heavy roller may compact a soil toa 
greater density than the light roller. For extremely light or heavy rollers, 
little compaction may be obtained. 

While it is true that great compacted density can be obtained only by the 
employment of heavy rollers, better efficiency of compaction can be attained 
by “stage compaction.” “Stage compaction” is the process of compacting 
first by the trips of light weight rollers, followed by the use of heavy rollers. 
The variation of efficiency and density of “stage compaction” with the number 
of passes are shown in solid lines in Fig. 2-A and Fig. 2-B respectively. If 
an exceedingly heavy roller is used in order to obtain high density, the louse 
soil is likely to be incapable of supporting the roller at the start. This diffi- 
culty can be overcome by “stage compaction” by using lighter roller first un- 
til the soil density is high enough to support the heavy roller. Two or even 
more stages can be used to advantage. 

In practice, if the job is of such size that two or more rollers may be used, 
it will be more economical and advantageous to use rollers which exert dif- 
ferent pressures instead of rollers which exert the same pressure. 

Furthermore, it is believed that a heavy roller can be designed so that its 
pressure for compaction may be increased gradually with the number of 
passes. With such a roller in being, the “stage compaction” method would 
not be necessary and the limit of compaction would be greatly extended. 


Moisture Content 


As is generally known, the influence of moisture content upon density is 
one of the most important factors in soil compaction. The moisture content 
affects principally the resistance to relative movement of the fines of soil. 

Each clay particle may be coated with a film of water molecules, one or 
more molecules in thickness. The first layer of adsorbed molecules, pos- 
sessing definite orientation, is ice solid and is exceedingly viscous. As more 
water becomes available, the adsorbed water layer becomes thicker and, at 
the same time, the degree of orientation decreases toward the outer layer un- 
til a thickness is reached in which orientation is nil. The viscosity of the 
water film decreases with the degree of orientation of water molecules. 
Further addition of water to the soil simply provides a fluid lubricant of low 
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viscosity and further separates the distance between clay particles. The 
frictional resistance to relative displacement of clay particles is hydrody- 
namic, as previously stated. 

The effect of adding water is to decrease viscosity and increase the dis- 
tance between particles. According to Equation 5, the frictional resistance 

to relative motion between particles is thus reduced. For the same compac- 
tive energy, therefore, the dry compacted density varies with the moisture 
content. 

The typical dry density versus moisture content curves are shown in Fig. 3. 
The validity of Equation 5 in explaining the effect of moisture content upon 
density is quite obvious. For cohesive or moderately cohesive soil, the com- 
pacted dry density is nearly directly proportional to the moisture content un- i 
til the “optimum” is reached. For materials containing very little or no fines, j 
the moisture content has almost no effect upon the density. The gravelly soil 
with just enough fines to fill the spaces between gravel particles is the most 
sensitive to the moisture content, while clay is the least sensitive. This is 
because the frictional resistance against movement varies directly as surface 
areas, assuming other factors are the same. A unit weight of clay possesses 
far more total surface area than an equal weight of gravelly soil. The “opti- 
mum” moisture content for cohesionless sand is saturation. This is not due 
to any “lubrication” action of water, but rather due to the reduction of effec- 
tive pressure between sand grains by hydrostatic pressure. 


Optimum Moisture Content and Optimum Compacting Pressure 


A soil mass is generally considered to consist of solid substances, and 
voids which are filled, or partly filled, with liquids and gases. For our study, 
it may be regarded as composed of solid particles, water and air. In com- 
paction of cohesive or slightly cohesive soil, assuming that soil particles and 
water are incompressible and that the water does not have time to be drained Te 
away, the air voids are the only remaining portion of the soil that is subjected 
to modification. The compaction of a soil mass is, therefore, the process of 7 
forcing soil particles closer together and thus of reducing its air voids. 
The variations of the relationships among solid soil particles, water, and 
air under the same compactive effort are shown by Fig. 4. The optimum : 
moisture content is clearly indicated in the diagram as the moisture content 
corresponding to the point of maximum solid substance. The soil mass which 
has a water content beyond the optimum will have a lower density because the 
excess water occupies the voids or space which would otherwise have been 
filled with solid particles. The air void approaches approximately a constant 
value after the optimum water content is reached. It represents the volume 
of air trapped in the soil during compaction. 
The term “optimum moisture content” is used rather loosely among engi- i 
neers.: It is often forgotten that optimum moisture content of a coil type is 
that for one compactive pressure or one compactive effort only. Each com- 
pactive pressure has its own optimum moisture content. 
There are in common use numerous types of rollers with different pres- 
sures. It is impossible to predetermine the optimum moisture content of a 
soil without knowing the type of the roller that is to be used in the field. The 
so-called Proctor and Modified A.A.S.H.O. laboratory compaction methods 
are valuable standards for use in studying the compaction behavior of a soil, 
but they should not be used indiscriminatively as means to determine the : 
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optimum moisture content of soil for field compaction. For good laboratory 
control of field compaction, it is often necessary to vary one or more of the 
details of the standard Proctor compacting method such as the weight of the 
hammer, the height and number of hammer blows, the thickness of layers and 
even the size of compacting mold, in order to obtain a more nearly exact re- 
lationship between optimum moisture content of the material and the roller 
to be used. 

There are often situations in which the natural existing moisture content of 
the soil to be compacted is such that its modification for efficient compaction 
from an adopted roller is economically impractical. For example, in attempt- 
ing to construct roller earth dams with exceedingly wet borrow material, it 
may be infeasible to dry the soil because cf the large volume of materials in- 
volved and the delay in construction that may occur. Since the moisture con- 
tent under such conditions is nearly a fixe:i value, it follows that the “optimum 
compacting pressure” must be determined. This is a trial and error pro- 
cedure in which the pressure of the roller is varied until its field optimum 
water content nearly coincides with the existing water content. 


Speed of Rollers 


Dynamic resistance of soil is a well recognized phenomena in pile driving. 
It may also exist to a certain degree in soil compaction during the rapid pene- 
tration of the tamping feet of a sheepsfoot rcller or of the highly inflated rub- 
ber tire of the pneumatic tire roller. If fast enough, the load will be resisted 
not only by the static friction and cohesion, but also by the viscosity of the 
soil. 

Dynamic resistance of soil is a complex problem and relatively little is 
known about the subject. However, the results of a number of experiments 
on the action of dynamic forces, especially in connection with the application 
of dynamic pile-driving formula, are valuable in arriving at qualitative under- 
standing of the problems. Figure 1-A indicates the relationship of penetration 
versus resistance and penetration versus velocity.(5) It may seem that the 
resistance will rise to a maximum value, which is principally dynamic re- 
sistance, at the beginning of the penetration and then drop rapidly to an al- 
most constant value, which is static resistance. Figure 1-B illustrates the 
variation of dynamic resistance with velocity. The area under the resistance- 
velocity curve represents the loss of energy. The shape of the curve is sig- 
nificant in that the area increases at a much more rapid rate than the in- 
crease in velocity. It is therefore more efficient to reduce the velocity of 
penetration if the energy to be consumed is constant. 

However, since the dynamic resistance is a direct function of velocity, and 
the energy input is also a direct function of velocity, equations may be written: 


Mx = [Rds 
Mvdv=Rds 


S=MJ§ vdv (° 
Je ftv) 6) 
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in which R = soil resistance 
S = settlement 
M and V are mass and velocity, respectively, as previously stated. 
In practice, the equation R = f (v) for the soil under consideration is too 
complex to be determined. Therefore, Equation 6 is of academic value only. 
Nevertheless, there may be formed a concept that there exists the most bene- 
ficial velocity to obtain maximum settlement, or the optimum velocity for 
compaction. 


Size of Rollers 


The influence of the size of rollers on compaction is similar to the in- 
fluence that the size of a footing has on settlement. This influence varies 
with the types of soil. For cohesionless soil, the difference in sizes of the 
loading area has little effect on compaction so long as the pressure remains 
the same. However, for the very heavy roller, shear failure will start sooner 
under small loading areas than under large loading areas. For clay, rollers 
with small contact areas are more efficient for compaction than those with 
large contact areas, in that there will be more volumetric strain for the 
small loaded area under the same pressure intensity, provided, of course, 
the soil is not loaded to failure. 

For soil types between cohesionless sand and cohesive clay, the pressure 
strain curve of a roller is as shown in Fig. 5, (7) assuming it is similar to 
that of a footing. It is seen that if the pressure is constant and does not ex- 
ceed q, the large area will have produced more strain than the small area. 
Therefore, the roller with the larger contact area will be more efficient than 
that with the smaller contact area. Actual tests in which sheepefoot rollers 
with varying sizes of tamping feet were used, have substantiatec the above 
statement.(8) Nevertheless, attention is called to the fact that by Feeping all 
the other things constant, the energy input for compaction is directly propor- 
tional to the area of the tamping feet. 


Thickness of Lift 


gl 


The thickness of the hard crust that develops on top of loose soil under 
concentrations of heavy traffic on uncompacted fill is a good indication of the 
depth limit of compaction of the particular traffic vehicles. As the soil depth 
increases, the pressure under load decreases and the bearing capacity of soil 
increases. The yield, and thus the compacted density, will decrease with 
depth. Based upon the conception of loading and resistance, the effective 
compacting depth will depend upon the size of the loaded width. The 50 ton 
heavy pneumatic tire roller as now manufactured and used, may have a tire 
width as large as 30 inches. Its effective compacting depth may well exceed 

6 inches, which is the arbitrary thickness of compacted layer adopted for 
many specifications of rolling regardless of type of roller. On the other hand, 
the fact that the tamping foot of sheepsfoot roller may be in the order of 3 
inches in diameter, is the reason for considering using a much thinner lift for 
rolling by sheepsfoot roller than by pneumatic tire rollers. It is believed 
that, due to lack of established relationship and concept between types of roll- 
er and thickness of lift, the practice of adopting any arbitrary thickness of 
soil layer leads to one of the major sources of wasted energy during com- 
paction. 
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Number of Roller Passes 


There are two primary considerations for determining the number of 
roller passes for compaction; the increasing in soil resistance with the num- 
ber of rolling trips, and the foot contact area of the roller. As already ex- 
plained, each succeeding roller pass is made on a soil having higher resist- 
ance to deformation than that under the previous rolling. For a chosen 
roller on a certain type of soil, there is a definite number of passes which 
will give the maximum return in increasing density for the energy spent. 
(See Fig. 2-A) Familiar among compaction engineers is the “diminishing 
return” in increasing density when the number of passes is great. However, 
one simple but fundamental concept should be recognized: by using a certain 
roller, there is a maximum density of soil that can be compacted under any 
practical condition. Greater density may be obtained by using a heavier 
roller. (See Fig. 2-B) Too many passes beyond that for maximum return is 
' another principal source of “wasted” energy during compaction. 

It must also be realized that the foot contact area of one pass of a sheeps- 
foot roller is only about 5 to 17 percent of the ground area covered by the 
roller depending of course, on the type of roller, while the tire contact area 
of one pass of pneumatic roller is 100 percent of the ground covered. This 
fact may explain why it requires fewer passes for a pneumatic tire roller 
than for a sheepsfoot roller to produce a given density. 


Influence of Gravel Content on Compaction 


The laboratory Proctor test used for field compaction control is generally 
made on a fraction of the total sample passing through No. 4 sieve. Thus it 
is not based on a representative sample of the actual material. It is neces- 
sary to understand the effect of including the gravel fraction (material re- 
tained in No. 4 sieve) on the compaction in order to evaluate the field com- 
paction results based on laboratory control tests. This evaluation is of great 
importance in actual practice. For example, actual experience with earth 
dam construction by the writer’s firm, has shown that it is imperative to be 
able to clearly correlate the degree of compaction achieved in the field with 
the compaction results achieved in the laboratory. 

The theoretical total density of soil containing gravel can be obtained by 
the formula(9) 


in which 7, = total dry density of soil and gravel; 
Ys = dry density of soil (minus No. 4 size fraction); 
Yg = dry density of gravel (plus No. 4 size fraccion); 
p = percentage of gravel content expressed as a decimal. 


The dashed lines on Figure 6 show the variation of the total theoretical 
density, and the theoretical weights of gravel and soil versus the percentage 
of gravel content. Equation 7 is based on the assumptions that the soil density 
remains the same for different gravel content and that the unit weight of 


| 
‘ 
% Ys 


862-12 SM1 January, 1956 


gravel is the weight without any voids among individual gravel grains. In 
other words, it is assumed that the gravels merely float and displace an equal 
volume of soil. 

The validity of the above assumption should be noted. The gravels in the 
soil, if grouped together tend to form voids which may not be filled, or partial- 
ly filled with soil fines. If the gravel content increases, the chance of form- 
ing such voids by grouping also increases. The distribution should follow a 
statistic law. Therefore, the deviation of the assumption from the actual con- 
dition increases with the gravel content. This has been borne out by actual 
compaction experience. 

Equation 7 is theoretically correct when the gravel content is zero. Its 
applicability decreases as the gravel content increases. When the material 
is 100 percent gravel, the unit weight, according to Equation 7, is that for the 
solid rock or close to 155 pounds per cubic foot, which is obviously incorrect. 
The unit weight of 100 percent gravel is close to 110 pounds per cubic foot. 

The actual and corrected density and weight curves are shown in Fig. 6 in 
solid lines. The actual total density curve is the curve which may be ob- 
tained by field in-place density test. It deviates more and more from the 
total theoretical density as the gravel content increases, and approaches the 
actual unit weight of gravel. From the actual total density, the curves of the 
actual weights of gravel and soil can be computed. At any gravel content, the 
degree of compacting of the soil compared with that free from gravel is ob- 
tained by dividing the actual weight by the theoretical weight. 

For example, referring to Figure 6, the dry density of soil fraction passing 
through No. 4 sieve is 114 pounds per cubic foot; the bulk specific gravity of 
gravel (plus No. 4 sieve sizes) is 2.45, from which the dry density of gravel 
without voids is found to be 2.45 x 62.4 or 153 pounds per cubic foot. At, say, 
40 percent gravel content, 


Actual total density = 122.5 lbs. per cu. ft. as determined by field in-place © 
density test; 
1 


Theoretical total density (from Eq. 7) = 0400.60 * 127 lbs. per cu. ft. 


153 114 
Theoretical weight of gravel = 127 x 0.4 = 50.8 lbs; 
Theoretical weight of soil = 127 (1-0.4) = 76.2 lbs; 
Actual weight of gravel = 122.5 x 0.4 = 49.0 lbs; 


Actual weight of soil = 122.5 (1-0.4) = 73.5 lbs; 
Degree of compaction of soil = 13.3 x 100 = 96.4% of theoretical weight 


Ratio of volume of voids to total volume of gravel, or porosity = 


Hence, it is seen that, at 40 percent gravel content the soil has been compact- 
ed to only 96.4 percent of the soil density free from gravel and the gravel has 
3.5 percent voids. 

By repeating the similar procedures for materials with different gravel 
content, curves as shown in Fig. 7, may be obtained. The curve on the density 
of soil shows the hindrance to the compaction of fine material by the presence 
of gravels. The curve on the porosity of soil indicates the volume of voids 
formed by clusters of gravels. It is interesting to note that they are directly 
proportional. 
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SUMMARY AND CONCLUSION 


The work done during compaction of a soil mass is either stored as inter- 
nal energy or converted into heat in overcoming the frictional resistance that 
arises from the relative movement of soil particles. The latter may result 
in the decrease of air voids and increase in soil density. The most efficient 
method of compaction is therefore that which will perform work to accom- 
plish the desired objective; that is, the readjustment of soil particles accom- 
panied by decreasing voids. 

Soil compaction is essentially the process of reaching equilibrium between 
dynamic loading and soil resistance. It depends upon the intensity of pressure 
of loading rather than the total loading. High soil density may be obtained by 
using “stage compaction” with the use of light pressure roller first, followed 
by heavy roller. The effect of moisture content is to reduce the frictional re- 
sistance to the relative movement of soil particles. It is only effective on fine 
soil particles. Optimum pressure should be determined by taking into account 
that the soil moisture content may be fixed. Large contact area of roller will 
give better compaction of moderately cohesive soil than roller with smaller 
contact area, assuming the contact pressures are the same. The thickness of 
lift for compaction should be based on the contact width of rollers. Fora 
certain type of roller on a given soil, there is a definite number of passes 
which will give the maximum return in increased density for energy spent. 
Finally, the influence of gravel content upon compaction is important because 
it hinders the compaction of the fine grained soil fractions. 
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ST. LOUIS DISTRICT 


C. I. Mansur,! M. ASCE, and R. I. Kaufman? 
(Proc. Paper 864) 


Levee System and Underseepage Problem 


Soil conditions in the alluvial valley of the Mississippi River are such that 
seepage beneath levees along the river has always posed a serious problera 
during high waters. In fact, some crevasses in levees along the Mississippi 
River have occurred as a result of sand boils and subsurface piping. 

Prior to 1938 many of the levees along the Mississippi River in the St. 
Louis District, Corps of Engineers, were below project flood grade, and only 
relatively low heads could develop against them. Although some underseepage 
had been experienced along these levees, little serious trouble had occurred, 
mainly because of their low height. By authority of the Flood Control Acts of 
1936 and 1938 the St. Louis District initiated a program to complete the levee 
systems and to raise the grade of existing levees so that they could withstand 
the project flood. This program is now approximately 90 per cent accom- 
plished. (A plan of the present levee systems along the Mississippi River 
from Alton to Gale, Illinois, is shown in fig. 1.) ‘The increased height of the 
levees has created a problem which did not exist previously, inasmuch as 
considerable head can now be developed against the levees, and their safety 
against underseepage has become questionable. 

Although it was realized that underseepage would be a problem along the 
new levee system, the control of underseepage was not specifically provided 
for in the design of the new levees or enlargements. After experiences in the 
1951 and 1952 floods in the Omaha and Kansas City Districts of the Corps of 
Engineers, the Office, Chief of Engineers, and the St. Louis District requested 
the Waterways Experiment Station to make a comprehensive investigation of 
the seepage problem along the new levee system and to design the necessary 
control measures. This investigation was initiated in October 1952 and 
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essentially completed in June 1955. Approximately 70 per cent of the control 
measures indicated by these studies have been completed and considerable 
additional work is under contract. 


Development of Underseepage and Sand Boils 


Whenever a levee is subjected to a differential hydrostatic head of water 
as a result of river stages being higher than the adjacent land, seepage enters 
the pervious substratum beneath the levee through the bed of the river, river- 
side berrow pits, and/or the riverside top stratum, and creates an artesian 
head and hydraulic gradient in the sand stratum under the levee. This gradi- 
ent causes a flow of seepage beneath and landward of the levee as illustrated 
in fig. 2. The seepage from the pervious substratum that emerges at or land- 
ward of the levee toe is generally termed underseepage; concentrated under- 
seepage that carries sand up to the surface through a more or less oper chan- 
nel in the top stratum is called a sand boil. Piping is the active erosion of 
sand or other soil from below the ground surface as a result of substratum 
pressure and concentration of seepage in localized channels. An example of 
some small sand boils which occurred during the 1937 flood is shown in fig. 3. 

If the hydrostatic pressure in the pervious substratum landward of a levee 
becomes greater than the submerged weight of the top stratum, the excess 
pressure will cause heaving of the top blanket, or will cause it to rupture at 
one or more weak spots with a resulting concentration of seepage flow in the 
form of sand boils. Where the foundation and top strata are heterogeneous in 
character, as is usually the case, seepage tends to appear at localized spots 
instead of causing the entire top stratum to heave or become “quick.” 

The hydraulic gradient required to cause heaving is called “the critical 
hydraulic gradient;” it is the ratio of the submerged unit weight of the soil to 
the unit weight of water and is usually about 0.85. Any tendency for the hy- 
draulic gradient to increase above the “critical” gradient only causes addi- 
tional sand boils or increased percolation. In nature, high exit gradients and 
concentrations of seepage are usually found along the landside toe of the levee, 
at thin or weak spots in the top stratum, and adjacent to clay-filled swales or 
channels. 

Where seepage is concentrated to the extent that turbulent flow is created, 
the flow will cause erosion in the top stratum and development of a channel 
down into the underlying silts and fine sands which frequently exist immedi- 
ately beneath the base of the top stratum. As the channel increases in size 
and/or length there is a progressively greater concentration of flow into it, 
with a still greater tendency for erosion to progress beneath the levee. This 
is especially true if the top stratum is cohesive and the underlying soils are 
susceptible to erosion. Shrinkage cracks, root holes, and holes made by man 
or burrowing animals form ready channels in the top stratum for development 
of localized flows and attendant piping. Not only is underseepage a hazard 
from the standpoint of underground erosion, but it may also saturate and re- 
duce the stability of the landside slope of a levee to such an extent that the 
slope may slough with attendant dangers to the levee. 


Geology of Alluvial Valley and Its Relation to Underseepage 


Detailed geological studies of several sites along levees in the Lower 
Mississippi Valley where underseepage had been a problem during previous 


| 
7 
1 


January, 1956 


q 


FIG. 2. CROSS SECTION OF LEVEE AND ALLUVIAL VALLEY 


FIG. 3. SAND BOILS IN SACK SUBLEVEE BASIN NEAR FRIARS POINT, MISS. 
1937 HIGHWATER 
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high waters show a definite correlation between the distribution of alluvial 
deposits of sand, silt, and clay and the location and occurrence of underseep- 
age and sand boils. Thus, before it was possible to make an intelligent in- 
vestigation of underseepage and to design control measures, an understanding 
of the alluvial fill in the Upper Mississippi Valley from Alton to Gale, Illinois, 
was necessary. The geological studies made for this project included field 
reconnaissances, studies of aerial mosaics, and a review of logs of borings. 

A few miles below St. Louis, the Mississippi River cuts across an eastward 
extension of the Ozark Plateaus and for a distance of about 100 miles flows in 
a narrow (3 to 10 miles wide), rock-walled, alluviated valley incised into 
ancient, predominately limestone bedrock. The river emerges from the 
plateaus through a narrow gorge known as Thebes Gap into the broad Lower 
Mississippi Valley. The depth of the alluvial deposits in the St. Louis Dis~- 
trict is rather variable, ranging from about 75 to 200 ft with an average depth 
of about 125 ft. 

Generally, the river is adjacent to the bluffs on the Missouri side with the 
floodplain existing east of the river in the State of Illinois. Tributary streams 
originate in the bluffs on the east, flow through the floodplain, and empty into 
the Mississippi River at frequent intervals between Alton and Gale. The levee 
system, on the left bank, in this area consists of levees which originate at the 
bluff, run parallel to the tributary streams toward the Mississippi River, and 
then along the Mississippi River to the next tributary stream (see fig. 1). 
Thus the alluvial valley is subdivided into small areas, each with its own en- 
circling levees. 

Melting of the glaciers during late or postglacial times caused the en- 
trenched valley to become filled with a series of sandy gravels, sands, silts, 
and clays which can be grouped into two broad units: (a) a sand and gravel 
substratum, and (b) a fine-grained top stratum. As these units are of basic 
importance to the underseepage problem, they are described in some detail 
below. 


The Pervious Substratum 


During the gradual rise of sea level accompanying the final retreat of 
glacial ice, the Mississippi River became an overloaded braided stream in 
which large quantities of gravel-bearing sands were deposited. As sea level 
continued to rise and the deposits on the entrenched valley floor continued to 
thicken, stream slopes were progressively reduced and both the quantity and 
grain size of the material transported by the river decreased. The gravel- 
bearing sands were succeeded by coarse sands grading upward into progres- 
sively finer materials and terminating in deposits of very fine sand. This up- 
ward gradation from coarse to fine materials reflects the gradual adjustment 
between the transporting capacity of the river and the available load. Fine to 
very fine sands were not deposited until sea level had reached essentially its 
present stand and the river began to change from a braided to a meandering 
stream. The substratum sands are quite pervious and have a high seepage- 
carrying capacity as subsequently discussed. 


The Top Stratum 


The predominantly fine-grained top stratum deposits have been laid down 
by the Mississippi River as it meandered across the floodplain. Meander belt 
deposits are characterized by marked lateral and vertical discontinuities, and 
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by wide variation in grain size and permeability. They are therefore of 
primary impor:ance in the localization of underseepage and the design of 
control measures. Meander belt deposits may be grouped as point bar, chan- 
nel fi/l, natural levee, and backswamp deposits. Such deposits can usuaiiy be 
identified in the field and on aerial mosaics. 


Point Bar Deposits 

As meander loops increase in radius by erosion, deposition takes place on 
the inside bank where low sandy ridges are built up with intervening elongated 
depressions which usually become filled with fine-grained deposits. Ridges 
and swales formed in this manner are known as point bars. Long sandy ridges 
separated by clay- and silt-filled depressions more or less paralleling the 
margins of former river courses are also prominent features of shifting river 
valleys and are especially common in the St. Louis District. These features 
appear to have been formed in migrating channels rather than on points and 
might well be referred to as channel bars and slough fillings. The upper part 
of these ridge-and-swale deposits are usually covered with clay silts and silty 
clays which are laid down during gradual migration of the river channel from 
its former banks. The slough fillings normally are about 10 to 20 ft. deep. 


Channel-fill Deposits 

When the river abandons its former course as a result of a cutoff, the 
central and lower portions of the old cut-off channel usually become filled 
with silts and clays which are relatively impermeable. Such deposits in the 
area studied may be 40 to 80 ft. deep. 


Natural Levees 

When the river overtops its banks the water spreads out, the velocity is 
checked, and deposition of a portion of the sedimentary load results. In this 
manner, long ridges known as natural levees are formed on the outside of 
meander loops and along both banks in straight reaches. Natural levee de- 
posits found in the Upper Mississippi Valley appear to be mainly sandy in 
composition and are seldom over 200 yards in width and about 5 ft. high. 
Clearly distinct natural levees are relatively rare in the St. Louis District. 


Backswamp Deposits 

Low-lying areas on the landside of natural levees are known as back~ 
swamps. These areas receive only quiet flood waters and as a result the 
sediments deposited consist of fairly uniform silts and clays. Backswamp 
deposits create an almost impervious block to the emergence of subsurface 
seepage. There do not appear to be many, if any, true backswamp deposits 
in the St. Louis District; however, some thick clay top stratum deposits were 
found adjacent to the valley walls which resemble backswamp deposits. 


Effect of Alluvial Deposits on Underseepage 


The emergence of seepage landward of a levee is influenced by: (a) con- 
figuration of geological features, such as swale fillings and channel fillings 
and their alignment relative to the levee; (b) characteristics andsthickness of 
the top stratum; (c) cracks or fissures formed by drying and other natural 
causes; (d) borrow pits, post holes, seismic shot holes, ditches, and other 
works of man; and (e) decay of roots, uprooting of trees, animal burrows, 
crayfish holes, and other organic agencies. The severity of the underseepage 
along a reach of levee is frequently dependent upon the configuration of geo- 
logical features in that area. The greatest concentration of seepage always 
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occurs along the edges of swales and the landside levee toe as illustrated in 
fig. 4. 


Field Investigations 


An extensive field investigation was made to ascertain the need for under- 
seepage control measures and to obtain data on which to base their design. 
The field investigation consisted of a reconnaissance of all areas adjacent to 
the levees, study of aerial mosaics, an intensive boring program, a geophysi- 
cal investigation to determine the depth of the alluvial valley, and field pump- 
ing tests. 


Field Reconnaissance 


Field reconnaissances were made along the levees for the purpose of not- 
ing features appearing on the ground surface that might affect the underseep- 
age problem. The presence of geological features such as swales and sloughs 
and their orientation with respect to the levee were noted. The locations and 
extent of man-made features such as drainage ditches, stock ponds, and bor- 
row pits were ascertained. The layout of new borings, made in the office on 
the basis of available soils data, levee heights, and a study of aerial mosaics, 
was checked in the field to insure adequate exploration of those features which 
would affect the underseepage problem and the design of control measures. 
The condition and depths of riverside borrow pits also were noted.. 


Aerial Mosaics and Maps 


As no complete maps of the new levee system were available to whic’ the 
new borings, piezometers, and control measures could be referenced, a new 
fo‘io of photo-maps was prepared showing all Mississippi River levees and 
f'.ood-control works in the area. These new maps will be used during flood- 
« ontrol operations and to record underseepage observations during future 
:loods. A portion of ene of these maps, including the Harrisonville Levee 
District between stations 795 and 820, is shown in fig. 5. 


Boring Program 


At the beginning of the investigation the existing soils information was re- 
viewed and a boring program planned to delineate soil conditions more pre- 
cisely in potentially critical areas. Shallow borings to define the character 
and thickness of the top stratum along the landside toe of the levees were 
located on about 200- to 250-ft. centers and at critical locations. Some shal- 
low borings also were made farther landward of the levee where necessary to 
delineate pertinent geological features. A few borings were made in riverside 
borrow nits to determine the character and thickness of the remaining top 
stratum in these pits. In general, blanket borings were made with hand augers 
whenever water-table conditions would permit. Deeper borings (60 to 80 ft.) 
were made on 1000- to 1500-ft. centers to determine the characteristics of 
the subsurface sands for both the design and installation of relief wells in 
reaches where it was thought that control measures might be necessary. 

Some borings were extended to rock to determine the thickness of the pervious 
aquifer and to check geophysical determinations of the depth to rock. 

Most cf the deep borings were made with a 2-in. split spoon sampler using 
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drilling mud to keep the hole open. Samples were usually taken on about 5-ft. 
intervals and at major changes in soil strata. A few undisturbed samples of 
the sand foundation were obtained by use of 3-in. Shelby tubing and drilling 
mud, (4) The reason for obtaining the undisturbed samples was to compare 
grain size of the sands as determined from such samples with that of sam- 
ples obtained with split spoon and bailer samplers. Another reason for tak- 
ing a few undisturbed sand samples was to study the stratification of the sand. 
A typical soil profile developed from field investigations is shown for a short 
reach of levee in the Harrisonville Levee District in fig. 6. 


Geophysical Investigations 


Because of the difficulty and cost of ascertaining the depth to rock by 
means of borings, it was decided to obtain this information by means of 
geophysical explorations. 

At the beginning of the investigation the depths to rock were determined by 
the electrical resistivity method. In analyzing the data obtained it became 
apparent that the depth to rock determined by the resistivity method was ex- 
tremely sensitive to the method of interpretation and did not check depths 
determined by borings too well. As a result this method was abandoned in 
favor of the seismic apparatus. 

Initial checks were made with a seismograph at locations where the 
depth to rock had been determined from deep borings, and from excavation 
for the Chain of Rocks Canal Lock. Good agreement was obtained and, 
therefore, the seismic method was used for the remainder of the investiga- 
tion. 


Field Pumping Tests 


The design of relief wells along the levees in the St. Louis District was 
based on the assumption that the screens for the wells would penetrate the 
pervious aquifer underlying the levees to a depth which would result in per- 
formance of a well system comparable to that which would be obtained by 50 
per cent penetration of a homogeneous, isotropic aquifer. As the pervious 
strata underlying the levees along the Mississippi River generally become 
more pervious with depth, special pumping tests were performed on wells 
fully penetrating the pervious aquifer to determine how deep the screens for 
the wells should be installed to obtain an effective penetration of 50 per cent. 
The over-all permeability of the sand foundation at the cites of the test wells 
was also determined from these tests. The horizontal permeability of the 
various sand strata encountered at the test wells was determined by measur- 
ing the flow into the weli screen from each sand stratum as previously deline- 
ated by a boring. 

The field pumping tests consisted of pumping the wells at three different 
rates of flow, measuring the drawdown in the wells, and observing the draw- 
down in adjacent piezometers and wells. The flow in the well screen was 
measured at changes in strata by means of a calibrated well flow meter. (2) 
Head losses through the filter and well screen and inside the well were mea- 
sured by piezometers installed in the periphery of the gravel filter around the 
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well screen. The over-all horizontal permeability Ky of the sand aquifer and 
of the individual sand strata was computed from the formula for artesian 
radial flow to a single well. 

Grain-size curves of typical sands taken from a boring at one of the test 
wells (FC-105) together with the gradation of the filter gravel used around the 
well screen are shown in fig. 7. The permeabilities of the different sand 
strata encountered at well FC-105 as determined from laboratory tests and 
field pumping tests are shown in fig. 8. 

The flow of the test wells was essentially artesian because of the existence 
of upper and lower impervious strata bounding the principal aquifer, and the 
fact that the water in the well was not drawn down below the top of the main 
sand stratum. As would be expected with artesian flow, drawdown in the test 
wells stabilized rapidly at a constant rate of pumping; approximately 80 per 
cent stabilization was achieved in 15 to 30 minutes of pumping; full stabiliza- 
tion in 2 hours of pumping. Drawdown curves to the test wells as plotted on 
semilog graph paper resulted in straight lines, another indication of artesian 
flow to the test wells. The radius of influence of the different test wells 
ranged from about 500 to 1000 ft. The specific yield from the test wells 
ranged from about 150 to 375 gpm per ft. of drawdown at the well. Head loss 
through the filter and wooden well screen was quite small, amounting only to 
about 0.10 to 0.25 ft. for a flow through the well screen of 10 gpm per ft. of 
screen. In analyzing the well test data, it was necessary to take into account 
the hydraulic head losses in the well screen and riser pipe. 

The results of the field pumping tests may be summarized as follows: 
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WELL BORING 
WFC-I05 WB 
4 
WELL FC. 108 LOCATED AT LEVEE STA 301+25, FT. CHARTRES LEVEE 
TRICT, APPROX. 45 RIVER MILES SOUTH OF ST. LOUIS. 
2. PERMEABIL:; TIES SHOWN BY BAR GRAPH WERE COMPUTED FROM WELL METER 
© PIEZOMETER READINGS. 
. PERMEABILITIES SHOWN BY @ AND © WERE OBTAINED FROM LABORATORY 
ESTS ON BORING AND WELL DRILLING SAMPLES. RESPECTIVELY, AND WERE AD. 
WSTED TO THE ESTIMATED NATURAL VOID RATIO BY THE FORMULA K,,=K, ($2) 
NO CORRECTED TO T +20° C. 
. FIGURES TO LEFT OF BORINGS ARE Dio IN MM. OBTAINED FROM BORING 
SAMPLES: THOSE TO RIGHT OF BORING /ROM WE.L ORILLING SAMPLES. 
S. CIRCLED NUMBERS REFER TO GRAIN SIZE CURVES SHOWN ON FIG. 7. 
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FIG. 8. COEFFICIENT OF PERMEABILITY AND EFFECTIVE GRAIN SIZE OF 
INDIVIDUAL SAND STRATA - WELL FC-105 


a) In general, little agreement was found between permeabilities as deter- 
mined in the laboratory on remolded samples and those obtained from the 
field pumping tests. Ther? is no reason why the permeabilities should agree, 
because the aquifer is stratified and lenses of coarse sand and fine gravel 
exist. Generally, the field permeabilities for any given strata exceed by 2 to 
4 times the permeability as normally determined in the laboratory. 

b) A correlation obtained between effective grain size Dj9 vs ky as deter- 
mined from flow from individual sand strata is shown in fig. 9. The scatter 
of points can probably be attributed to variation in uniformity of the grain-size 
curves and the fact that the D,,,’s used in the plot are based on the average of 
relatively few D;9’s for each sand stratum tested. In the absence of pumping 
test data, the graph shown in fig. 9 may be of use to engineers concerned with 
seepage, water supply, or irrigation problems in the Mississippi River Valley. 

c) The permeability of the principal seepage-carrying stratum at the test 
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well sites ranged from 1500 to 3000 x 1074 cm per sec. The permeabilities 
of individual sand strata varied from 200 to 6000 x 10-4 cm per sec. Pump- 
ing tests on relief wells subsequently installed throughout the St. Louis Dis- 
trict indicated that the permeability of the principal seepage-carrying sand 
stratum generally ranges from about 1000 to 2000 x 10-4 cm per sec. 

d) In order to achieve an effective penetration of 50 per cent of the princi- 
pal water-carrying stratum the wel! screen should penetrate approximately 
60 per cent of the depth of the principal sand aquifer. 


Laboratory Tests 


Laboratory studies included visual classification of all samples, grain-size 
determinations and permeability tests on numerous representative samples, 
and filter tests on the gravel filter used around the well screens. 


Classification 

Samples were classified both visually and, when in doubt, by means of 
mechanical analyses. The classification used was based on the Unified Soil 
Classification System. 


Grain-size Determinations 

Mechanical analyses were made on any samples of the top stratum whose 
classification was questionable, and on numerous sand samples obtained from 
the foundation for the purpose of determining the D;g and Dgs sizes. These 
values were used in delineating the principal seepage-carrying stratum 
(sands with Djg > 0.15 mm) and determining those strata in which weli 
screens should not be set because of the possibility of the foundation sands 


migrating into and through the filter gravel. Although tests showed that the 
filter would satisfactorily drain sands with a Dgs < 0.15 mm, well screens, 
in general, were not set in sands with Dgs5 < 0.25 mm. 


Permeability Tests 

Permeability tests were conducted on a number of remolded samples of 
sand. The coefficients of permeability as determined in the laboratory were 
adjusted to a temperature of 20°C and to the estimated natural void ratio of 
the sands. The natural void ratio of the sand was estimated from the Dgo of 
the cand.(4) Coefficients of permeability as determined from laboratory tests 
on aplit spoon samples of the alluvial sands ranged from about 100 to 2000 x 
10°-* cm per sec. 

In connection with the permeability tests made on samples obtained with 
split spoon samplers in mudded holes it is pointed out that such samples fre- 
quently are contaminated by drilling mud. Therefore, it was necessary to 
wash all such samples thoroughly before performing any permeability tests. 
Some contamination of samples taken with a Shelby tube in a mudded hole was 
also noted. 


Filter Tests 

Several laboratory tests were also made to check the ability of the filter 
gravel specified for relief wells to protect foundation sands against piping in- 
to the well. The filter sample tested was on the coarse side of the gradation 
band specified. The thickness of filter tested was 6 in., the minimum required 
in the field. 

Three relatively uniform sands typical of the finer foundation sands en- 
countered in the St. Louis District having Dg, values of 0.20, 0.15, and 
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0.12 mm were used in the test. The tests proved that the filter would drain 
safely all of the samples tested at gradients less than 8. The maximum 
gradient existing at the interface between the sand foundation and filter dur- 
ing either pumping tests or flood periods when wells are in operation will 
probably not exceed 1. Therefore, it is concluded that if the filter material 

is satisfactorily placed and has a minimum thickness of 6 in. it should protect 
any sand with a Dgs coarser than 0.16 mm. 


Methods of Underseepage Control 


There are several methods which may be used for controlling seepage be- 
neath levees. The choice depends on a number of factors, including the 
character of the foundation, cost, permanency, availability of right-of-way, 
maintenance, and disposal of seepage water. Methods that may be used to 
control underseepage are landside berms, relief wells, impervious riverside 
blankets, sublevees, cutoffs, drainage blankets, and drainage trenches. Only 
the first two of the above methods were considered practicable for the control 
of underseepage along the levees in the St. Louis District. 


Landside Berms 

Landside berms have been used extensively for the control of underseepage. 
A berm serves to strengthen the top stratum, increase the base width of the 
levee, reinforce the landward slope of the levee against sloughing during high 
water, and to move the emergence of seepage or sand boils and piping out 
from the levee proper. Landside berms were used along levees in the St. 
Louis project wherever they were the most economical and where the top 
stratum landward of a levee was very thin or nonexistent and relief wells 
would have had to be installed excessively close together. Landside berms 
were also used where a levee was relatively low but yet it appeared that a 
certain amount of protection against underseepage and possible sloughing of 
the landside slope was desirable. 


Relief Wells 

One method of controlling underseepage is to tap the underlying pervious 
strata with a series of relief wells which will provide pressure relief and 
controlled seepage outlets that offer little resistance to flow and, at the same 
time, will prevent erosion of the foundation soils. This method has been used 
successfully by the Corps of Engineers at a number of dams and along numer- 
ous levees underlain by a pervious foundation. Relief wells offer an advantage 
as compared to gravel toes, pervious blankets, and other surface measures 
where the foundation consists of stratified deposits of pervious materials, in 
that they penetrate the more pervious strata in which pressure relief is 
necessary. 

With proper spacing and penetration, relief wells will reduce hydrostatic 
pressures landward of levees underlain by impervious foundations for a wide 
range of seepage entrances, foundation stratification, and landward top strata. 
In addition, seepage which formerly emerged landward of a levee in an uncon- 
trolled manner without relief wells will be materially reduced, although the 
total of well flow plus seepage with wells will be increased somewhat. (5,6) 


Impervious Riverside Blanket 

Frequently the natural top stratum riverside of levees has been removed 
as a result of borrow operations for construction of the levee, thereby expos- 
ing the pervious substratum and providing a ready source of seepage entry. 


| 
| 

| 

| 
| 
| 
| 
Ned 
| 
| 
| 

ks 
Ay 

| 
i 
| 
7 
a 


Paper 864 MANSUR - KAUFMAN 17 


No artificial impervious blankets were placed riverward of the levees in the 
St. Louis District as a means of seepage control. However, where the river 
had scoured deep holes or channels in the riverside borrow pits these holes 
or channels have been filled, either by dredged or hauled material, and abatis 
dikes have been constructed to promote silting and the growth of willows in 
the pits (see fig. 20, page 50). 


Sublevees 

Areas known to be susceptible to dangerous underseepage may be encircled 
by one or more sublevees in which water can be impounded during high water 
to reduce the net effective head acting in this critical area. Construction of 
sublevees requires considerable right-of-way landward of a levee, skillful 
operation during high water to insure proper balancing of heads to prevent 
dangerous sand boils or blowouts, adequate spillways and adequate section 
to insure that the sublevee itself will not fail during maximum flood stages 
owing to a head of water in the sublevee basin. For these reasons they were 
not considered practicable as permanent underseepage control measures in 
the St. Louis District. 


Cutoffs 

Where practicable, the most positive method of underseepage control is to 
cut off all pervious strata beneath a levee by means of an impervious barrier 
which will eliminate both excess substratum pressures and the problem of 
seepage water landward of a levee. However, completely cutting off pervious 
strata 80 to 200 ft. deep along extensive reaches of levees is not economically 
feasible. The installation of partially penetrating cutoffs will not reduce seep- 
age and excess pressures significantly unless the cutoff penetrates 95 per 
cent or more of the pervious aquifer. This fact has been demonstrated by 
both mathematical and theoretical studies.(7) No deep cutoffs were attempted 
along the levees in the St. Louis District for these reasons. However, a few 
shallow cutoffs along the riverside toe of the levee have been constructed to 
cut off relativeiy thin layers of either natural levee or crevasse sands which 
lie immediately under the base of the levee and which in turn are underlain 
by more impervious strata. 


Drainage Blankets 

Drainage blankets may be used for controlling underseepage where the 
levee is built on exposed sands and gravels of fairly homogeneous character. 
However, they are not effective for controlling seepage in deep substrata, 
where the drain would be underlain by impervious top strata, or where strati- 
fied fine sands would exist between the drain and the deeper, more pervious 
sand. No drainage blankets were used to control deep underseepage in the 
St. Louis District. 


Drainage Trenches 

Drainage trenches may be used to control underseepage where the pervious 
foundation is of limited depth so that the trench will substantially penetrate 
the formation. Where the pervious foundation is deep a drainage trench of 
any practicable depth would attract only a small portion of underseepage, its 
effect would be local, and detrimental underseepage would bypass the trench. 
Because of the depth of the pervious substratum along the Mississippi River 
levees, drainage trenches were not considered feasible for this project. 
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Seepage control measures were considered necessary wherever the esti- 
mated head beneath the top stratum at the landside toe of the levee or existing 
berm, at flood stages equal to the net grade of the levee, would create an ex- 
cessively high upward gradient through the top stratum. In agricultural areas 
where it appeared that the upward gradient, i, would become equal to or ex- 
ceed 0.85, seepage control measures were considered necessary. Where the 
gradient was between 0.67 and 0.85, piezometers were installed to record the 
substratum pressure developing during future high-water periods; from these 
data decisions will be made later as to the need for control measures. Where 
the gradient was less than 0.67, no control measures or piezometers were 
considered necessary. In industrial areas, control measures were considered 
necessary wherever i was expected to exceed 0.67, and were designed so 
that the gradient at the levee toe would not exceed 0.50. 

The head, h_, beneath the top stratum at the landside toe of existing levees 
was computed from the following formula: 


wherein 


H = the design net head on the levee 
s = the distance from the landside toe of the levee to the effective 
seepage source 
Xq = the distance from the landside toe to the effective seepage exit 


The design net head for riverfront levees was taken as the difference in ele- 
vation between net grade of the levee and the average ground surface land- 
ward of the levee. Along the flank levees the design flood stage was taken as 
the elevation of the net grade of the riverfront levee at the intersection of the 
flank and riverfront levees. The distance to the effective source of seepage, 
S, was estimated from field reconnaissances of borrow pits and the river, 
boring data, and piezometric data obtained at similar sites along levees in 
the Lower Mississippi River Valley. For riverfront levees, values of s 
ranging from 600 to 1000 ft. were used for analytical and design purposes. 
Along flank levees where the seepage entrance is restricted to the bed of the 
tributary stream, values of s ranging from 800 to 1500 ft. were used. The 
distance from the landside levee toe to the effective seepage exit, x3, was 
computed from the following formula by Bennett,(8) assuming the top stratum 
level and infinite in landward extent. 


where 


c = a constant 
k,/ k, = retio of horizontal permeability of the foundation to vertical perme- 
ability of the top stratum 


—— 


— 
x 
4 
ian 
} 
i 
4 
| | 
i 
Age 
| 
| 
7 
| 
| 
| 
| 


Paper 864 MANSUR - KAUFMAN 19 


z = top stratum thickness in feet 
d = effective thickness of aquifer in feet 


To simplify the preparation of design charts it was assumed that d = 80 ft; 
variations in d between 50 and 120 ft. do not greatly affect the values of hp 
or Xq. Evaluation of the top stratum thickness and the permeability ratio was 
predlanted largely upon the thickness of the clayey portion of the top stratum. 
The value of z was obtained by transforming the thickness of each layer 
making up the top stratum into an equivalent thickness of material of constant 
vertical permeability. Where the top stratum was predominantly silty with 
less than 5 ft. of clay, kg/kp) was assumed equal to 400. Where the top stratum 
consisted of 5 ft. or more of fat clay, ke /ky, was assumed equal to 1000. The 
above values of k,/k, were based on analyses of piezometric data from simi- 
lar sites in the Lower Mississippi River Valley. By substituting xg as given 
by equation (2) in equation (1), the head h, can be expressed in terms of H 
and z for a given s, ke / , and d. Design curves used for computation of 
hydrostatic head at the landside toe of the levee without relief wells were 
developed from equations (1) and (2); a typical set of curves for s = 800 ft. is 
shown on fig. 10. 

Methods and typical graphs used for designing relief well systems and 
seepage berms are outlined in the following paragraphs. 


Relief Wells 


Relief well systems were designed so that the gradient midway between 
wells would not exceed 0.67 in agricultural areas and 0.50 in industrial areas. 
The effective well radius, ry, was taken as 0.9 ft.; i.e., the distance to the 
periphery of the gravel filter. Design curves were developed from Barron’s 
equations(9) for an infinite line of fully penetrating wells with a semipervious 
top stratum by computing the head midway between wells and then multiplying 
this head by the ratio of the head computed for 50 per cent and fully penetrat- 
ing wells as obtained from formulas by Jervis and Muskat,{10) respectively, 
for an impervious top stratum. As the computed head midway between wells 
does not include hydraulic head losses in the well, they were computed from 
the estimated flow from the well system and from hydraulic formulas, well 
tank tests, and observations made on relief wells at Trotters, Mississippi, (6) 
and added to the computed head midway between wells. The well spacing, ao, 
was then adjusted so that the final head between wells, including hydraulic 
head losses, would not exceed the design values stated above. Design curves 
were developed for kp /Kiy = 400 and various sources of seepage. Where k/Kpy 
= 1000, a, = 0.9 (a, for ke/k, = 400). A typical set of design curves showing 
the relationship between the computed well spacing (a,), adjusted for the head 
losses in the wells, and top stratum thickness for various heads on a levee is 
shown in fig. 11. 

The flow from relief well systems was obtained from Barron’s formulas 
for fully penetrating wells multiplied by the ratio of flow from 50 per cent 
penetrating wells to that for 100 per cent penetrating wells for an impervious 
top stratum as determined by Jervis and Muskat, respectively. The relation- 
ship between well flow, Q,/H, for various top stratum thicknesses and well 
spacings is shown in fig. 12 for ke = 1000 x 1074 cm per sec and s = 800 ft. 
From results of pumping tests conducted to date it now appears that the 
average ky equals about 1200 to 1500 x 10-4 cm per sec for the pervious sub- 
stratum in the Upper Mississippi River Valley. 
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In a short, finite line of wells the heads midway between wells exceed those 
obtained for an infinite line of relief wells both at the center and near the ends 
of the well system. As numerous well systems were fairly short (less than 
1500 ft. in length) it was necessary to reduce the well spacing computed for 
an infinite line of wells so that heads midway between wells in a short line of 
wells would not create an i in excess of the design value. The ratio of the 
head midway between wells at the center of finite weil systems to the head 
between wells in an infinite line of wells is shown on fig. 13 for various well 
spacings and exit lengths. Where a short line of wells was required, the spac- 
ing required was computed for an infinite line of wells, and this spacing then 
was divided by the ratio indicated in fig. 13. In any finite line of wells of con- 
stant penetration and spacing, the head midway between wells near the ends of 
the system exceeds that at the center of the system. Thus, at the ends of both 
short and long well systems the wells were generally made deeper to provide 
additional penetration of the pervious substratum so as to obtain the same 
head reduction as in the central part of the well line. 

The average well spacing for a given reach of levee was determined in the 
office from design curves similar to those shown on figs. 10-13, taking cog- 
nizance of the effect of geological features on seepage conditions. However, 
the location of each well was checked in the field and adjusted where necessary 
so that the wells would fit natural topographic features. 


Seepage Berms 


Seepage berms were designed by an empirical method developed by the 
Kansas City District, CE, based on electrical analogy model studies in which 
it was assumed that the permeability of the berm would be the same as that 
of the top stratum. Adding a berm on the landside of a levee restricts the 
natural relief of pressure as a result of natural seepage through the top stra- 
tum and thus increases the hydrostatic head at the levee toe with respect to 
the original ground surface. This increase in head depends upon the length, 
thickness, and permeability of the berm. The electrical analogy model studies 
indicated that the effect of the berm on landside substratum pressures would 
be about the same as though the impervious base width of the levee had been 
increased a distance equivalent to 50 per cent of the length of the berm. The 
nomenclature used in the design of seepage berms is given in fig. 14. The 
procedure used in designing seepage berms in the St. Louis District is out- 
lined below. 

a) The head, h_, at the landside toe of the levee without a berm was deter- 
mined from equation (1). (The allowable i, with a berm was taken as 0.60 in 
agricultural levee districts and 0.50 in industrial districts.) 

b) The allowable head at the berm toe, h,, was taken as 0.85 z in agricul- 
tural levee districts and 0.67 z in industrial districts. 

c) The first trial length of berm, xj, was computed as that length of berm 
where the head, h,, at the berm toe would be equal to h,; this length was de- 
termined by means of Bennett’s (8) blanket formulas. 

d) The value of x;/2 was then added to s , and the head, No, at the point 


x,/2 distance from the landside toe of the levee was computed from formula 
(1), which results in the following equation: 
Hx 
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HYDROSTATIC HEAD AT LANDSIDE TOE (NO WELLS) IN FEET (h,) 


Note: Relief wells steve om 


THICKNESS OF TOP STRATUM IN FEET (2) 


FIG. 10. COMPUTED HYDROSTATIC HEADS AT LANDSIDE TOE OF LEVEES - 
NO RELIEF WELLS 
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THICKNESS 2F TOP STRATUM IN FEET (2) 


FIG. 11. COMPUTED WELL SPACING ADJUSTED FOR HEAD LOSS IN WELL 
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the Performance of Relief wells," Proceedings — 
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LENGTH OF WELL REACH-FEET 


FIG, 13. RATIO OF HEAD MIDWAY BETWEEN WELLS AT CENTER OF A FINITE WELL 
SYSTEM TO HEAD MIDWAY BETWEEN WELLS IN AN INFINITE LINE OF WELLS 
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e) A new distance, x41, from the assumed impermeable section (x;/2 dis- 
tance from the landside toe) was computed by means of blanket formulas so 
that the head h, at distance xj, would then be equal to the allowable head h,. 

f) The required berm length X is then x1/2 + x34. To simplify the above 
computations it was found expedient to represent the blanket formulas graphi- 
cally as shown in fig. 15, which shows the relationship between the ratio of 
head, hy, at x distance landward of the levee to the head, h,, at the toe and the 
distance x divided by the exit length, x3. By selecting h, and substituting 
this value for h, and knowing h, and xg, the distance x at which h, = h, can 
readily be found from the above figures. 

To compute the required thickness of the berm it was necessary to esti- 
mate the head, hd, at the actual levee toe with the berm in place. This was 
accomplished by means of the following formula: 


x 
+ x,) 

hese 
tx, 


The average upward gradient through top stratum and berm was expressed by 
the following equation: 


h-t 
z+t 


i= 


where t = berm thickness. From this equation the required berm thickness 
t at the levee toe can be found as follows: 


° 


where i, is the allowable upward gradient at the levee toe. 

The design of a seepage berm required for the reach of levee between sta- 
tions 810 and 817 in the Harrisonville Levee District, had relief wells not 
been installed, is illustrated in fig. 14. 

The authors have developed an equation for determining the required length 
X of a berm from Bennett’s blanket formulas and unpublished equations by 


DISTANCE FROM CENTER LINE OF LEVEE In FEET 
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FIG. 14. NOMENCLATURE FOR DESIGN OF SEEPAGE BERMS 
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FIG. 15. RATIO BETWEEN HEAD LANDWARD OF LEVEE AND HEAD AT LANDSIDE 
TOE OF LEVEE 


Bennett for design of seepage berms. This equation is as follows: 


A 24 (2 +r) 


wherein 
A= 6+3sc (r+1) 
r=i,/iy 
ig = allowable upward gradient at landside toe of the levee, ; 


ij = allowable upward gradient at landside toe of the berm, 


and H, hg, s , and c are as described previously in this paper. It should be 
noted that two values of X are obtained from equation (4). As a result of the 
sign convention adopted in developing the sore, the negative value of X 
gives the required berm length. The head, h+ , at the toe of the levee with 
the berm in place can be computed from the Fades equation: 


phen, fr-cx+(2$*) (5) 


The required thickness of the berm at the levee toe is given by equation (3a). 
The above equations developed by the authors indicate the following required 
dimensions for the seepage berm illustrated in fig. 14: X = 175 ft., t = 2.5 ft. 
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Installation and Construction of Seepage Control Measures 


Relief Wells 


The relief wells along the levees in the St. Louis District have been and 
are being installed by the reverse rotary method. This procedure for drilling 
holes for wells in sand is basically a suction dredging method in which the 
material in the hole is removed by a suction pipe; the walls of the hole are 
supported by seepage forces acting against a thin film of fine-grained soil on 
the walls of the hole created by maintaining a head of water in the hole several 
feet above the water table. A film of fine-grained soil is deposited on the 
walls of the hole by the drilling water as water and soil from the suction pipe 
are circulated through a sump pit in which the sand settles out and from which 
the water containing the fine-grained particles flows back into the hole. Suc- 
cessful drilling by the reverse rotary method requires that the water table be 
about 7 ft. or more below the ground surface and that there be an adequate 
supply of drilling water. 

A typical reverse rotary drilling unit is normally equipped with a large- 
capacity centrifugal pump, 6-in.-diameter drill pipe, and a bit similar in ap- 
pearance to the cutterhead of a dredge. Where cobbles larger than about 
3-3/4 in. are encountered it is necessary to provide this type of equipment 
with a rock trap. One of the reverse rotary drill rigs used on the St. Louis 
project is shown in fig. 16. Holes for relief wells also may be advanced by 
the reverse rotary method using compressed air or a jet eductor system, 
neither of which requires the use of rock traps. 

A typical relief well consists of a screen section, riser pipe, gravel filter, 
sand backfill from the top of gravel filter to an elevation 10 ft. below finished 
ground surface, and concrete backfill from the top of the sand backfill to 
ground surface. The riser and screen are of 8-in. I.D. wood pipe. The screen 
is perforated with 3/16-in. slots + 1/32 in. wide. The bottom of the well 
screen is closed with a wood plug. Tops of relief wells discharging on the 
ground surface are protected by metal well guards. A schematic drawing of 
a typical relief well is shown in fig. 17. 

As the boring for a relief well is advanced, samples of the foundation soil 
are obtained at 2-ft. intervals by catching samples of the effluent from the 
drill rig in a bucket. The purpose of this sampling is to determine the depth 
at which the screen section of the well can be started and also to locate strata 
of silts, silty sands, and very fine sands through which unslotted sections of 
pipe are used instead of slotted screen. During or immediately after drilling 
operations the well screen and riser pipe are assembled and guides attached 
to the pipe to keep it centered in the hole as the filter gravel is placed. (The 
hole for each well is overdrilled by at least 4 ft. so as to provide a space at 
the bottom for filter gravel which may become segregated when the tremie 
pipe is first filled.) The filter gravel is placed only after the screen and 
riser have been lowered to the correct depth and the top of the riser pipe set 
about 4 in. above the natural ground surface. The filter gravel is placed by 
first lowering a perforated tremie to the bottom of the hole and filling it with 
filter gravel. The tremie then is slowly raised to allow the filter gravel to 
run out of the bottom as it is fed into the top. The tremie pipe is kept filled 
at all times with the filter gravel being kept above the water surface so as to 
prevent any segregation of the gravel as it is placed. The filter gravel is 
placed to at least 4 ft. above the top of the screen section so as to insure the 
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FIG. 16. WELL DRILLING RIG IN OPERATION 
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existence of filter above the top of the screen after the well has been surged 
and developed. 

After the filter is placed, material which has entered the well during the 
placing operations is removed and the well is surged and pumped to remove 
drilling mud or other fines from the filter gravel. Surging and development 
pumping are started within 2 hours after the filter has been placed and are 
continued until the amount of material pulled through the screen between surg- 
ing operations is less than 0.2 ft. deep in the well. Surging is usually accom- 
plished by lowering and raising a surge block made of a heavy rubber disk 
between two steel disks mounted on a steel rod. The rubber disk has a diame- 
ter approximately 1 in. less than the inside diameter of the well and the steel 
disks are 1 in. smaller than the rubber disk. The surge block is raised and 
lowered at a rate of travel of approximately 5 ft. per sec. 

Each well is subjected to a pumping test after completion of surging so as 
to achieve a drawdown of 5 ft. in the well or a flow of at least 500 gpm per 
minute. The inflow of sand into a well during the pumping test is carefully 
checked for each well to prove its stability. This is accomplished by pumping 
the flow into a carefully baffled 1000-gallon sedimentation tank for 15 minutes. 
The sand in the tank is collected and measured, and the change in amount of 
sand in the bottom of the well ascertained during this interval. The pumping 
test is continued until the rate of inflow of sand into the well is less than 1 
pint per hour. A relief well installation report is shown in fig. 18. Fig. 6 
illustrates a completed reach of well system, showing the soil profile, spacing 
of wells, head on levee, location of well screen and filter gravel, and includes 
a bar graph of the specific yields of each well. A photograph of a completed 
well system is shown in fig. 19. 

The top of each well is protected against backflow of muddy surface water 
by means of a specially designed rubber gasket and simple check valve as 
shown in fig. 17. Accelerated full-scale laboratory tests have indicated this 
check valve to be most effective in preventing backflow of muddy water into 
the well under adverse conditions. 

During periods of relatively low stages on the levees, relief wells will flow 
somewhat more than natural seepage. In agricultural areas each well is pro- 
vided with a plastic sleeve which will raise the discharge elevation of the well 
to either 1 or 1.5 ft. above the natural ground surface. These low sleeves or 
standpipes will prevent well flow at low river heads on the levee, at which 
stages no pressure relief is necessary. As soon as artesian pressure de- 
velops to the extent that water begins to spill over the top of the standpipes, 
they will be removed and the well system allowed to operate as originally de- 
signed. The flow from relief wells in highly developed areas will be pumped 
over the levee during high water. 

A systematic program of maintenance of relief wells and piezometers in- 
stalled in the St. Louis District is planned to insure their continued and effec- 
tive operation. 


January, 1956 


Piezometers 


Piezometers consisting of a 2-ft. section of 1-1/4-in.-diameter brass well- 
point, plastic coupling, and 1-1/4-in. galvanized iron riser pipe have been in- 
stalled in the pervious substratum immediately beneath the top stratum in al- 
most all well reaches to check the efficacy of the well systems during high 
water. Generally these piezometers are installed betweer. the wells near the 
end and in the middle of any sizeable reach of wells. At certain typical 


45 


a 
| 
| 
‘ 
4 
r 
: 
5 
~ 
Sid: 
i... 
: 
“4 
a 
| 
| 


Paper 864 MANSUR - KAUFMAN 29 
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DRILLING 


DATE AND /7 June ‘55 18 June 
4:30PM tle: 5:00AM 


wmoo: Leverve-Rofory From : 


tye: Wood stave, machine bonded , waite tam: 8-in. Pertortions: x 


Wood shove pipe, 
mre: Grove/ 


 Golvonized corrugeled me/a/ TYPE OF 
wei cum: pipe with wire vereen man: Concress 

Lrechmann 
cometon: LuAr Bros. Inc. contact orematon : Deer 
GROUND ELEY «Fine To 
AT WELL : 378.2 W:- Sa: 276 fine Sond: 27. 
TOP TOP OF 
Yum: 378.60 380.0 375.35 
TYPE OF 
wn: A/uminum plote 


FIMAL WELL INSTALLATION DATA 


«0 
Screen 
~ 
Cloy 
Design depth of wel/= 84 Cobbles encountered 


of depth of 79.3 Fh. Unable to obsain design depth. 


FIG. 18. RELIEF WELL INSTALLATION REPORT 


TO SAD 
ELEV TOP OF 
“QESMN” SCREEN OL WELL SCREEN 547.5 Riser 
fey TOP 
TOTAL SCREEN LENGTH : 39.1 OF PUTER: 368.8 
wen: 727 300.95 
DEPTH ELEY 
378.7 
Wey Smith 
| 
| 


SM 1 January, 1956 


locations in each levee district a few piezometers have been installed on a 
line perpendicular to the levee and beneath it for the purpose of estimating 
the effective source of seepage entry at those points and as a check on design 
assumptions and performance of the well systems. Piezometers were also 
installed in areas where the estimated upward gradient at the design flood was 
between 0.67 and 0.85. The screen for the piezometers have either 40-mesh 
brass screen or No. 18 slots. The tops of the piezometers are protected with 
a column of concrete 6 in. in diameter which extends from just below the 
piezometer cap to 30 in. below the ground surface. After each piezometer is 
installed it is pumped and tested to check its performance. The locations of 
some piezometers in plan are shown in fig. 5. 


Seepage Berms 


Seepage berms may be constructed by either hydraulic fill methods or by 
hauling. Berms constructed of sand are preferable but in general those con- 
structed to date in the St. Louis District are built of random soil. However, 
in no event should a berm be constructed of soil that, when placed, is more 
impervious than the natural top stratum. 
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JOURNAL 
SOIL MECHANICS AND FOUNDATIONS DIVISION 
Proceedings of the American Society of Civil Engineers 


STRUCTURAL RIGIDITY IN CALCULATING SETTLEMENTS 


Samuel Chamecki,! M. ASCE 
(Proc. Paper 865) 


1, INTRODUCTION 


The assumption, currently used in the calculation of the settlements, that 
the supports of statically indeterminate structures possess infinite facility to 
adapt themselves to the deformations of the foundation soil, leads to conclu- 
sions diverting from reality, especially in the case of high building structures. 

The continuity among its various members gives an actual structure con- 
siderable rigidity, with the result that the differential settlements become 
much lower than when calculated under the assumption above. 

The attempts to consider the rigidity of the structure in calculations of its 
foundation settlements, using processes based on assumptions differing from 
the real elastic behavior of the structure and neglecting the real rheological 
constants of the soil, yields results as uncertain as those derived from the 
simple application of empirical correction coefficients. 

The rigorous and general solution given herein to the problem of inter- 
action between structure and foundation soil is consistent in all respects with 
what is orthodox in both structural and foundation engineering. 

This simplified process results from the introduction of some simplifica- 
tions in the exact solution which, because of being evidently reasonable, need 
no further justification. 

The tasks of the structural and soil engineers have been kept separated so 
as not to divert them from their routine work. The structural engineer de- 
velops a column load chart with indication of coefficients of load transference, 
which are elastic constants of the entire structure; the soils engineer, in 
possession of these data, computes the foundation settlements, without use of 
any notion strange to his specialized field. 

As we will see hereunder, the relations established by means of the coef- 
ficients of load transference, make possible an exact consideration of the 
rigidity of the structure in calculating foundation settlements, without risk of 
being affected by modifications which are likely to occur in consequence of 
progress in Soil Mechanics. 


Note: Discussion open until May 1, 1956. Paper 865 is part of the copyrighted Journal 
of the Soil Mechanics and Foundations Division of the American Society of Civil En- 
gineers, Vol. 82 No. SM 1, January, 1956. 


1. Prof. of Structural Eng., University of Parana, Pasana, Brazil. 
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2. General Solution 


To illustrate the exposition by a most simple and definite example, we will 
consider the case of a two-span continuous beam, resting on a foundation soil 
containing a compressive clay layer, which will be responsible for almost the 
whole settlement (Fig. 1). 


2.1 Calculation of the settlements not considering the structural rigidity: 

Representing by Roj the reaction at any support i, without consideration of 
the unevenness of the supports (but with consideration of continuity), 8; the 
contact area of any footing plate i , the contact pressure between plate and 


: 


soil will be: 


and, in the case of our beam (Fig. 2a): 


Poc* 


Be Nj, a number which, multiplied by the contact pressure Poj under any 
footing j , yields the compression stress resulting at a depth y, at the ver- 
tical passing through the center of the footing i. The depth y corresponds 
to the horizontal plane where the stresses represent the average one from 
which results the total deformation of the layer. 

Nii is a function of the area, shape and location of footing j , of the depth 
y and of the rheological characteristics of the soil. 

The total stress oi (at the vertical of the center of footing i, at the 
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depth y) will be: 
he *Mi PoitZ Poy | [ies]  ..-(2) 


In the case of spread foundations, the number Nj; corresponds to the num- 
ber of divisions, multiplied by the respective values of Newmark’s graphic 
(Fig. 2b), occupied by footing j; the center of the graphic corresponding to 
the center of the footing i and all scale reiations being respected. In the 
case of our beam, we have (Fig. ic): 


Naa Pos +Nea Poc 
Jop=Nas PoatNes Pos +Nea Poc 
Toc = Nac Poa +Nac Pos +Nec Poc 


The settlement at the vertical of any footing i, will be given by the ex- 
pression: 


T- being the modulus of deformability of the soil for the known value of stress 


In the case of normally consolidated clay, the settlement will result frorn 
the “void-ratio - stress-log.” curve (Fig. 2d), where: 


In the case of our beam, we have (Fig. 2e): 


H, 
Naa Poat Nga Poe +Nea Poc) 


H, 
4og= rome Nas Poa + Nae Poe + Nes Poc) 


ral Nac +Neac Nec Poc) 


2.2 Influence of structural rigidity: The settlements QO, (Fig. 2e), calcu- 


lated under the assumption of the footings not being mutually connected and 
therefore capable infinitely to adapt themselves to the deformation of the soil, 
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cannot correspond to reality. The rigidity of the structure modifies the pat- a 
tern of the settlements, generally resulting in lower differential settlements. ¢ 
Be the real settlements of our beam, the ones illustrated on Fig. 3a. They me 


will provoke in the structure resistances to differential settlements which will 
modify the reactions Rp at the supports. 


: 
8 
a 
F19. 3 
a In order to determine the influence of the settlements on the reactions at 
the supports, we will superpose the effects, considering separately each sup- - 
port’s settlement. 
If support A undergoes a settlement of unit value, it will create in A, B 
and C, the respective vertical reactions Qaa, Qap and Qac (Fig. 3b), the 
sum of which must be zero to satisfy the condition of static equilibrium: 
Thus, -Q = 2 Qij; i.e., the reaction at the footing i , undergoing a unitary vot 
settlement, is equal and opposite to the sum of reactions created at the other te 
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footings. In other words: support i, when settling a unit «quantity, gets re- m 
lieved of Q;; Which is transferred to the other supports. 

Therefore we will name the QQ, coefficients of load transference between 
the supports, the first index indicates the procedure, i.e., the footing settling 
a unit value, and the second index, the local, i.e., the support to which reac- ' 
tion Q corresponds, or which receives the transferred load. 

Qij is read: Qfrom i to j. 

From the Maxwell’s reciprocity theorem, we know that Qj = Qi - 

As we will see later, it will be by means of the coefficients of load trans- 
ference that we completely will characterize the influence of the rigidity of 
the structure in the calculation of the foundation settlements. 

If the support A undergoes a total settlement Aa instead of a unit once, 
the resulting reactions in A, B and C will be, respectively (Fig. 3c): 


In the same way, in figures 3d and 3e, the reactions created by the separate 
settlements 2) g and va\ c appear. These reactions are always given in the 
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An 


form of a product of a coefficient of load transference and a corresponding ; 
settlement (index of settlement being the same as the first index of Q). 7 


Thus, if Roa is the reaction at the support A (Fig. 2a), without considera- 
tion of the unevenness of the supports, and if -%, ASA; + Opa 4 g and 
-Qc4 4, are the respective influences on that reaction of the settlements 
and Ac (Fig. 3c,d,e), we have in our beam, considering the 
(Fig. 4a): 


settlements 


and, in the same way: 


Rs =Rogt Gas Os Agr 


Re= Roc A, A; Ac 


Therefore, considering the rigidity of the structure, we now have the reac- 
tions at the supports given in function of the unknown settlements. 

The reacti~ns Ro are known and the coefficients of load transference Q 
are elastic constants of the entire structure which can be easily calculated 
by the processes of statically indeterminate structures (and correspond to 
the influence lines of reactions at the supports, for unit settlements).2 

For a support i, we will have: 


Ri = Roi Vii t2 ..- (5) 


2. Samuel Chamecki - “The influence lines in plane and three-dimensional 
continuous frameworks” ~- Editora Cientifica. Av. Erasmo Braga, 299 - 
8° andar. Rio de Janeiro, Brazil. 
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2.3 Calculation of the settlements considering the structural rigidity: We 


will obtain hereafter the set of simultaneous equations making possible the 
calculation of the unknown settlements, with consideration of the rigidity of 
the structure. 

If the influence of the rigidity of the structure on the settlements is im- 
plicit in relations of the form (5), expressing the reactions at the supports 
in function of the unknown settlements, we will repeat the operations de- 
scribed in 2.1, to ubtain the contact pressure p; and the vertical compressive 
stresses G; , in order to obtain finally, the wanted set of equations. 

Using expression (5) for Rj at a generical support i, we will have for the 
contact pressure at the support: 


fig. 4 


For our beam, we have (Fig. 4a): 


Gas Gee Ae dc 
Ps 


~ oe Gac 4a Ac 
Pe 
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> Re (6) 
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For the compressive stress at a support i, we will have: 


Applied to our beam (Fig. 4b), we have: 
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Ta = NaaPa +Naa Pe +Nea Pe ; 
05 = Nas Pa +Nas Pe tNes Pe 
= Nac Par Nec Pe + Nec Pe 
The settlement at a support i will result from: 
where the modulus of deformability Cj has to be expressed in function of | 
G; , for it varies with it, because Hooke’s law is not valid for soils, and the 
value of Or is not known, as in the case of the expression (3). Thus: 


Function C =f, (7) can be obtained for any soil and represents the 
derivative of the “stress - strain” curve, with respect to 7 (Fig. 5a,b). 

In our case of a normally consolidated clay, we have C = ae - av (Fig. 
5c); we can write down the set of the three equations with three unknown 
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quantities which, solved, yields the settlements at the supports Ay, ’ As 
and Be of the beam (Fig. 4c), calculating considering its rigidity. 
The equations are, in case of our beam: 


oe Yon Re Sc 


A- + Nec Pe Nec Pc) 


ae & 


We have chosen the verticals passing through the centers of the footings 7 
as the supports of the vectors representing the unknown settlements. This 
procedure is justified by the following reasons: = 
I) we avoid the necessity of inter or extrapolation, disadvantages which 4 
appear when we choose points not coinciding with the centers of the footings. es) 


Thus we may consider the special conditions encountered under each footing, 
i.e., thickness of underlying layers responsible for the settlements, modulus 
of deformability of soil, conditions of consolidation influenced by adjoining 
buildings, etc.; 

II) we accentuate the predominance of the principal diagonal of the coeffi- 
cient matrix of unknowns in the set of equations (8), which makes the solution 
possible by process of iteration, which, as we will show in paragraph 4, will 
be substituted by a sequence of simple operations, thus eliminating the neces- 
sity of establishing the equations. 

Being the NN and QQ with equal indexes generally of larger value than 
those of different indexes, it appears that the above mentioned diagonal is 
predominant. In fact, in the first equation, the unknown quantity 4, only is 
multiplied by Qa ,4Naaq, in the second one, only A g is multiplied by QanNpp 
and so on. 


3. Simplifications on the Calculation of the QQ 


The structural engineer has various ways of applying simplifications in 
the calculation of the coefficients of load transference, which depend on the 
importance of the structure and the precision required. 

In the case of building structures (Fig. 6a), we can calculate the QQ of a 
generic column i (Fig. 6b), considering the frame of an average story as 
consisting of a continuous beam passing through i, the columns of one story 
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above the one beneath, and the perpendicular beams on both sides of the con- e 
tinuous beam (Fig. 6c), and the other frame passing through i in a perpen- 
dicular direction, consisting of the same elements (Fig. 6d). The result is 
multiplied by the number of stories. 
Another simplification which can be done, is to consider the coefficients 
of load transference from each column to the adjoining ones. In this case, 
thanks to the reciprocity of the QQ, it will be sufficient to work out the cal- 
culation with the only consideration of the unit settlements of the columns of 
even numbers (in both directions). 


a) AW 
AY 
A 
Z 
d) 
fig.6 


In this condition, and in the case of our beam (Fig. 1), it would be suffi- 
cient to consider the support B as having settled a unit quantity and calculate 


the coefficients Qpa and 
The other coefficients therefore become known, because: 


Yes =-(Qsa 

Pas =-Qas 

he = Vac 

Ce 


and 
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4. Sequence of Calculation of the Settlements 


Let us consider briefly the substance of the method we propose to solve 
the problem of consideration of the rigidity of the structure in the calculation 
of the foundation settlements. 

In paragraph 2.2 we introduced the expression (5), which supplies the re- 
actions at the supports of the structure in function of the unknown settlements; 
in this expression appear the coefficients of load transference which are elas- 
tic constants of the structure and which determine the interaction between 
structure and foundation soil. 

In paragraph 2.3 we introduced the general expression (8), which leads to 
the set of simultaneous equations from the solution of which results the set- 
tlements, considering the effect of structural rigidity. 

The solution of this set of simultaneous equations being extremely labori- 
ous, we skirt around the difficulty by calculating the settlements at the ver- 
ticals of the centers of the columns which enables us (see 2.3) to consider the 
conditions of the subsoil particular to each point, and to solve the set by the 
process of iteration, the diagonal being predominant. 

As we will see hereafter, this process will be substituted by a sequence 
of operations, without necessity of establishing the equations properly, and 
making possible the final result by successive approximations. 

In order to keep separated the tasks of the structural and the soil engi- 
neer, we give hereunder briefly the process of the calculation of the settle- 
ments, with consideration of the rigidity of the structure, applicable to any 
type of soil and any type of structure. 


4.1 Task of the structural engineer: 

I) Calculation of the loads Ry of the columns (vertical reactions at the sup- 
ports), without considering unevenness of the supports. 

I) Calculation of the coefficients of load transference Q. 

Ill) Elaboration of the column location plan with indications of the respec- 
tive loads and coefficients Q, for supplying to the soils and foundation engi- 
neers. Using the simplifications of paragraph 3, this plan looks like figure 7, 
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because we considered only the coefficients Q for adjoining columns and be- 
cause Qij = Qii , by reciprocity. 


4.2 Task of the soil and foundation engineer: 

I) Design of foundation. 

II) Calculation of the settlements Ao , without considerjng the rigidity of 
the structure, i.e., for the loads Rg from paragraph 4.1 - I. 

Il) First correction due to the rigidity of the structure: 

a) calculation of the new loads values on the columns (expression 5), using 

the formerly obtained settlements J, (§ 4.2 - I); 

b) recalculation of the settlements, corresponding to the new column loads, 

according to identical procedure as in § 4.2 - II. These settlements re- 

sulting from a first rigidity correction, are represented by 4 r 

IV) successive corrections due to rigidity of the structure: calculation of 
the new loads on the columns by means of the settlements resulting from the 
prior correction and recalculation of the settlements so many times that they 
result to be equal in two successive approximations or with a difference equal 
or less than the desired approximation. These values will be the final settle- 
ments. 

The values of the settlements in the successive approximations, that cor- 
respond to the successive cycles of the resolution, by means of iteration, of 
the set of simultaneous equations (8), oscillate around and get nearer the final 
value, as can be seen in figure 8 (values Gp , A, P 47: etc.); in order to 
eliminate the necessity of a large number of corrections, one should apply a 
second correction, after the first one, by calculating the new loads of the 
columns for the average settlements (4, + 4, ) é (Fig. 8). 

In most practical cases there won’t be any necessity of further corrections, 
as we will see by means of an example we will give later on. 


5. Example 


To calculate the settlements of a structure represented in figure 9, with 
consideration of its rigidity. 

I) The column load values, without consideration of the unevenness, are 
given in figure 9c. 


3. The reaction of column D, in function of the unknown settlement will be: 
Ro =Roo + 


The coefficients of load transference Qnr = Qgp, for instance, multi- 
plied by the difference of settlements between the columns E and D (sub- 
tracting the smaller from the larger) supplies the load to be subtracted 
from that of the column which settles most and to be added to that which 
settles less. 

4. These settlements will be calculated by the usual processes of Soil Me- 
chanics, i.e., considering whether existing underlying layers particularly 
responsible for the settlements or being necessary to integrate the effect 
over the total layer; whether the soil deformation extends above the pile 
ends or under the footings, whether the soil is sand or clay, normally con- 
solidated or not. 

Our method appears to be of general application. 
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II) Calculation of the coefficients of load transference: using the simplifi- 
cation indicated in §3, perfectly assumable in the present case, it will be 
sufficient to calculate the Q coefficients for the frames passing through DEF 
and BEH, only for 4 = /, and multiplying the resulting values by 3 (3 
stories). 

III) In figure 10 the column load chart is given, together with the coeffi- 
cients of load transference which will be used to calculate the settlements 
with consideration of the rigidity of the structure. 

IV) Design of the foundations: the geotechnical survey indicates the exist- 
ence of an underlying uniform clay layer, which will be responsible for almost 
all of the settlement (Fig. 11a). 

The “void-ratio - log. of compression-stress” curve, of an average soil 
sample (Fig. 11c), shows the presence of a normally consolidated clay; in 
fact, for the point marking the preconsolidation, the stress J;= /0 ton/m? 
(1.024 ton/sq. ft.) corresponds to the load of the 6.00 m (19 ft. 8 in.) over- 
lying soil (Fig. 11b). 
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Fig. 


Making no considerations about the foundation’s design, which is of no 
concern here, we suppose that a foundation has been designed on isolated 
footings for a uniform contact load of Do = 40 ton/m? (4.10 ton/sq.ft.) (Fig. 
12), without consideration of the rigidity of the structure. 

V) Calculation of the settlements, without consideration of the rigidity of 
the structure: by means of Newmark’s graphic, we will calculate the vertical 
compression stresses in the center lines of the footings, in the average hori- 
zontal plane of the clay layer (elevation: -6.00 m or 19 ft. 8 in.). With the 
center of the graphic (in the suitable scale) in A, we have (the same results 
for C): 
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Fig. l 


May + Neg = 2037 ; = 0.008, 0.037; 
Meg =Q.013, Nog 20.004 ond =0.003. 
In B: Nag + Neg= 2010; 0.068; 
Neg = 0.074 Nog = 0.002 5 = 0.006 . 
D (equal toF, by symmetry): Nap tNep=2018; Ngp=O.006 ; 
Noo + = 071; Nep = 0.017; Neg # Nop =Q.0143 Nyp=0.006. 
In E : Nye t+ Neg =O0.006 Nye 5 Noe = ; 
Neg = 0.140 Neg + Nig = 0.006 3 Nyg = 0.070 
In_G (eguol to by symmetry): Nag = 9-003 ; 
Nog + Neg = 2024 3 Neg =O.0103 Nog 2050; Nyg=OOl2 . 
In_H: Naw Ney = 0.002 ; Nyy=O.004 Noy +Ney =0-010 ; 
Ney = 0.050 3 Now + Nyy =0.010 Nyy = 0.092 

The compressive stresses, using expression (2), will be: 


Goa = = (0.037 + 0.008 +Q037+0.013 +0.004+0,003)40=35.64 


The other stresses, calculated in a similar manner, are given in figure 13. 
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4068+/0 


With the data of the “ € - fog. 0” graph (Fig. 11c), we calculate the a 
corresponding settlements by using the expressions (3) and (4) and remem- 
bering that the settlements result from the final stresses (/- which equal 
the sum of the preconsolidation GT and of the q -stresses, due to the 
footings (Fig. 11b). 

Thus: 


E,-€ 09-€E 


for = 3.64+10=13.64 C/m?, we have: thus: 
= '.9-2.1 x 0.8785 =0.055=5.5 cm (2 . 


In a similar manner, we calculate the other settlements, which are given 
in figure 13b, from which the “curves of equal settlements” result (without 
consideration of the rigidity of the structure). 

IV) First correction due to the rigidity of the structure: by means of the 
settlements of figure 13b and the coefficients of load transference given in 
figure 10, the new column loads are as follows: 


Ria = Roa- Van Ma + + Ppa Ao = 
= 72-(2 670 + 801)0.06 + 801 x0.068+ 2670x0075 = 129 t 
or simply: 
Ria = Goa 44) + = 
= 72+ G01x0.028+2670 x 0.013 = 129 t . 
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Fig. 14 
The other reactions, the values of which are given in figure 14a, are cal- 4 
culated in a similar manner. ra 


Next, we calculate the contact pressures (expr. 6): 3 
: 
similar manner, the other pressures which are given in figure 14b. * 


The stresses /; are computed (expr. 7) by means of the same coefficients 
of Newmark’s graph used in the calculation above. 


Thus: 
= = 0.037 x 72 + 0.008 x49 + 0.037 x33 + 0.013 x 16 + 
+0.004 x64 +0003 x40= 5.11 t/m? . 
The results are given in figure 14c. 4 
The stresses being known, we calculate the settlements by means of the 2 
data of the graph of figure lic: iE 
21E 
for J=5.//+10 =15./1 t/m*?,wehave: €=0O87 , thus: 
= 19-2.1x 0.87 =0.073m 
The other settlements A ,? calculated in a similar manner, are given in ca 
figure 14d, 
VII) We know the settlements 4 o (without consideration of the rigidity of id 


the structure) and 4), (first correction due to the rigidity). 


We will make a second correction due to the rigidity of the structure, by Ro 
calculating the new column loads Ry, not for the above calculated settlements, Pe re 
but for the mean values in 


Ap (iz. 8) 
2 
In figure 15 are given the settlements representing mean values of the i 
settlements of figures 13b and 14d. ae 


qj 
i 
| 


865-18 January, 1956 


Ain (cm) 


Fig. 15 


The new column loads Roy the resulting contact pressure py, the stresses 
oz and the settlements 4 g , calculated in a way similar to the above opera- 


tions, are given in figure 16a, b, c, d. 


S59 4510 627+ G5. 
Pa (¢/m+) G+& 4, (cm) 


Fig. 16 


Comparing the settlements a Va with 4 m , we see there is no necessity 
of elaborating further corrections, because closer approximations do not 
make sense in Soil Mechanics. Consequently, the settlements 4, are the 
final calculated values, with consideraiion of the rigidity of the structure. 

In figure 16 the “curves of equal settlements” are drawn. 


== 


6. CONCLUSIONS 


There is no doubt that this resolution of the problem of the consideration 
of the rigidity of the structure in the calculation of its foundation settlements, 
is general, rigorous and simple. 

I) It is general, because it is applicable to any type of frame-structure, 
with any type of foundation, and resting on any type of foundation soil. 
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Moreover, we can consider the influence of the rigidity of the structure 
on the other components of the displacements of its foundation, i.e., not only 
the vertical one. It will be sufficient to calculate the corresponding compo- 
nents of the reactions (moment and horizontal reaction transference coeffi- 
cients). 

II) It is rigorous, because we have kept within what is orthodox in the two 
engineering fields of structures and soil mechanics. 

In order to calculate the coefficients of load transference, we may con- 
sider the structure in the three-dimensional space and the displacibility of 
its joints. 

The optional simplifications, suggested in §3, are evident to be reasonable. 

In all stages of the corrections, the settlements (i.e., soil deformations) 
are calculated by the usual processes of Soil Mechanics. In each stage it is 
sufficient to calculate the new column loads, by means of relation (5), which 
defines rigorously the interaction between structure and soil. 

This relation holds, even in the case of evolution and modification of the 
assumptions and processes for calculating soil deformations, for the various 
types of foundations. i 

III) It is simple because it modifies only to a small extent the calculation 
of the settlements in the case of no consideration of the rigidity of the struc- 
ture and because of practically not diverting the structural and soil engineers 
from their routines. 

Easy, too, is the analysis of the settlements and of their effects on the 
structure, ior various ages of it, considering the fact that we are working 
with elements which gives the influence of unit settlements. 
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SYNOPSIS 


The relative density of cohesionless soils and the bearing capacity of 
foundations in such soils can frequently be determined most conveniently 
from the results of penetration tests made on the site. From a comparison 
between standard (dynamic) and static penetration resistances on a number 
of sites, a simple approximate correiation has been obtained between both 
types of tests to estimate the relative density and angle of internal friction. 

The proposed relationships are applied to determine the ultimate bearing 
capacity of footings and the point resistance and skin friction of piles, and 
the results are compared with some field loading tests on plates and piles. 
The procedure can also be used to estimate the allowable load on spread and 
piled foundations in cohesionle.s soils under various conditions in practice, 


INTRODUCTION 


The ultimate bearing capacity of cohesive soils can generally be estimated 
from bearing capacity theory and the shearing strength of undisturbed sam- 
ples. For soils with little or no cohesion, however, the difficulty and expense 
of obtaining undisturbed samples makes estimates of the bearing capacity 
from penetration tests on the site frequently the most economical procedure. 
For this purpose dynamic or static penetration tests may be used. 

The most widely used penetration tests in U.S.A. and Canada are the stand- 
ard (dynamic) penetration tests in which a 2 in. outside diameter sampling 
spoon of 1 3/8 in. internal diameter is driven under an energy of 350 ft.-lb. 
into the ground ai the bottom of a borehole and the number of blows per foot 
of penetration is recorded.(1) On the other hand, in Europe cone penetration 
tests are the most common method and static tests are usually preferred in 
which a 60° cone of 1.4 in. base diameter is pushed into the ground at a slow 


Note: Discussion open until May 1, 1956. Paper 866 is part of the copyrighted Journal 
of the Soil Mechanics and Foundations Division of the American Society of Civil En- 
gineers, Vol. 82 No. SM 1, January, 1956. 

1. Head, Dept. of Civ. Eng., Nova Scotia Technical College, Halifax, Nova 
Scotia, Canada. 
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and fairly constant rate; both the cone penetration resistance and the skin 
friction on the casing tube surrounding the sounding rod are measured during 
the test.(2) In dynamic cone penetration tests the cone is driven into the 
ground under a constant energy and the number of blows per foot of penetra- 
tion is recorded; to avoid interference by skin friction the sounding rod is 
frequently encased. 

While static penetration methods are to be preferred, they suffer from 
the disadvantage that in dense and very dense soils a substantial reaction to 
jacking has to be provided or even the load capacity of customary equipment 
may be reached. Dynamic methods are cheaper and simpler to carry out and 
are fairly reliable in cohesionless soils. Cone penetration tests give a con- 
tinuous record and require additional boreholes for identification of the soil. 
Standard penetraticn tests, on the other hand, furnish disturbed samples for 
index property tests, but are often difficult to carry out satisfactorily under 
artesian ground water conditions. All penetration tests become unreliable 
as the maximum particle size approaches the diameter of the penetrometer 
or sampling spoon, and the lowest penetration resistance should be used if 
the maximum size of the soil particles exceeds about 1/2 in. 

Not only the procedure, but also the application of penetration tests differs 
considerably. Thus the results of standard penetration tests have been cor- 
related with the relative density of sands and the corresponding bearing ca- 
pacity of spread foundations. (3) The results of cone penetration tests, how- 
ever, have been correlated with the bearing capacity of piles. Static cone 
penetration tests can be extrapolated directly,(4) although for small penetra- 
tions of piles into a cohesionless stratum extrapolation by means of bearing 
capacity theory is to be preferred; (5) dynamic tests are usually extrapolated 
by means of a pile driving formula. (6) 

In order to make the penetration tests more widely applicable, it is there- 
fore proposed below to correlate the results of dynamic (in particular, stand- 
ard) penetration tests with the results of static penetration tests in cohesion- 
less soils so that the standard tests cannot only be used for estimates of the 
relative density and bearing capacity of spread foundations as previously, 
but also for estimates of the bearing capacity of piles, and so that the static 
tests cannot only be used for estimates of the bearing capacity of piles as 
previously, but also for estimates of the relative density and bearing capacity 
of spread foundations. 


Comparison Between Dynamic and Static Penetration Resistance 


Apart from some indirect comparisons between standard and static pene- 
tration tests on the basis of plate loading tests in shafts(7) or boreholes, (8) 
the only published accounts of direct comparisons relate to three sites in 
u.S.A.(9) and South America; (4) in all cases the variation of penetration re- 
sistance with depth was found to be similar for the standard and static pro- 
cedures. Further evidence of this similarity of resistance is provided by the 
results of tests at two sites in Eastern Canada. On each site tests were made 
at two locations to determine the variation of the standard penetration resist- 
ance, dynamic cone resistance, static cone resistance and static shaft friction 
(Figs. 1 and 2). The dynamic cone resistance (using a 60° cone of 2 in. base 
diameter on an uncased 1 5/8 in. diameter sounding rod driven under 350 ft.- 
lb. energy) is about twice the standard penetration resistance under the same 
energy, which may at least partly be explained by the area of the cone base 
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being about twice the net sectional area of the sampling spoon. With greater 
penetrations the skin friction on the sounding rod rapidly affects the mea- 
sured dynamic cone resistance and may lead to refusal in dense and very 
dense soils so that a casing tube should be provided. Apart from this limita- 
tion the dynamic and static cone resistances and.the standard penetration re- 
sistance vary in a similar manner with depth. The unit static skin friction, 
deduced from the measured total shaft friction, also follows a similar varia- 
tion with depth as the cone and spoon resistances. 

The results of the above-mentioned standard and static penetration tests 
have been plotted in Fig. 3, which shows that a reasonably good linear cor- 
relation exists between the standard penetration resistance and the static 
cone resistance. On the average 


Wc 4N 
where 4. static cone resistance (tons per sq. ft.) 


and N standard penetration resistance (number of blows per ft. pene- 
tration) 


equivalent penetration resistance (blows per ft.) 
15 +1/2(N - 15) for N>15 


in saturated very fine or silty sands. 

On the basis of the previous approximate relationship which has been sug- 
gested between the standard penetration resistance and the relative density 
of sands,(3) the above correlation between the static and standard penetration 
resistances may be used to suggest a similar approximate relationship be- 
tween the static cone resistance and the relative density of cohesionless soils, 
as given in Table 1. 

The present relationship between static cone resistance and relative 
density is similar to the somewhat more conservative tentative relationship 
deduced previously from a preliminary analysis of published data.(10) A 
close correlation between relative density and penetration resistance cannot 
be expected since the latter depends also on the soil type, overburden pres- 
sure, ground water conditions and, in the case of the standard test, also on 
the weight of the rods.(11) For the same _eason only an approximate relation- 
ship can be obtained between penetration resistance and the shearing strength 
or angle of internal friction on which this penetration resistance depends. 
Thus for different relative densities approximate angles of internal friction 
of cohesionless soils based on earlier suggestions (10,12,13) are given in the 
last column of Table 1. The lower values of the angles are safe limits for 
uniform clean sands and should be reduced by up to 5° for silty sands in the 
absence of shearing tests. The upper values are safe limits for well-graded 
sands and may be increased by up to 5° for sand-gravel mixtures if shearing 
tests are made on the material as a check. 


Bearing Capacity and Settlement of Spread Foundations 


The simplest type of spread foundation is a footing whose ultimate bearing 
capacity in homogeneous cohesionless soils may be expressed by 
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af = 7 + Y DN, 
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width of footing, 
depth of footing, 
unit weight of soil 


and Ny and N, are bearing capacity factors for a shallow footing and 
depend mainly onthe angle of internal friction of the soil.(5,14) Although the 
factors Ny and Ng depend on the assumptions made in different theories, (15) 
a safe expression is obtained by using 


<= 


G = + D/B) 
where q, = surface bearing capacity 
B 


From the relationship between N and @¢ given in Table 1, it can then be 
shown from Eq. 1 that, approximately, for a central vertical load on a footing 
on dry or moist sands 


Go = NB/10 (tons per sq. ft.) a a 


or = (tons per sq. ft.) 
where B is in feet, 


N = average standard (or equivalent) penetration resistance (blows 
per ft.) within depth B below base level, 


and Gc = average cone resistance (tons per sq. ft.) within that depth. 
Hence after substituting Eq. 4 into 3, 


= NB(1 + D/B)/10 (tons per sq. ft.) 
GcB (1 + D/B)/40 (tons per sq. ft.) time 
With a factor of safety of 3, the safe bearing pressure is then from Eq. 5 


or 


ds = NB(1 + D/B)/30 (tons per sq. ft.) eee 
or = q-B (1 + D/B)/120 (tons per sq. ft.) ee 


For silty sands the bearing capacities given in Eqs. 4 to 6 should be re- 
duced by up to one-half, while for sand-gravel mixtures the values may be 
increased up to twice in accordance with the above mentioned variation of N 
and @. Full submergence of cohesionless soils reduces the effective unit 
weight and, thus, the bearing capacities given in Eqs. 4 to 6 by about one- 
half. The bearing capacity is, however, not affected by a water table at a 
depth greater than about 1.5B below base level so that the bearing capacity 
yy jeterne positions of the water table can be interpolated according- 
ly. 

As a check of the proposed relationships between ultimate bearing capacity 
of footings and the penetration resistance of cohesionless soils, the results 
of plate loading tests(7,16,17) are given in Table 2 together with the corre- 
sponding estimates using Eqs. 4 and 5. It is found that the observed ultimate 
bearing capacities are conservative and are about twice the estimated values 
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in the case of 1 ft. wide footings. However, with increasing footing width the 
angle of internal friction at the ultimate bearing capacity decreases, (18) as 
also shown by the results of the loading tests in London (Table 2), so that the 
proposed relationships are considered reasonable for the footing widths com- 
monly used in practice. 

The allowable bearing pressure may be less than the safe bearing pressure 
if the settlement under the latter is excessive. It has been suggested(19) that 
a differential settlement of 3/4 in. can be tolerated by most ordinary struc- 
tures and that this movement will not be exceeded if the allowable bearing 
pressure causes 1 in. total settlement. From the chart for dry and moist 
sands,(3) it is found that this allowable bearing pressure is approximately 
given by 


N/8 (tons per sq. ft.) for BS4 ft. 
and qq = N(1 + 1/B)2/12 (tons per sq. ft.) for B>4ft. ...... (7b) 
or Ga = N/10, approximately, irrespective of B. 
Substituting Eq. 1 into 7a to 7c 

da = (tons per sq. ft.) for ft. 
and da = Gc(1 + 1/B)2/50 (tons per sq. ft.) for B>4 ft. RT 
or da = q,/40, approximately, irrespective of B. 


Comparison of Eqs. 6 and 7 shows that, irrespective of the relative density 
of the soil, the bearing pressure is governed by settlement considerations, 
i.e. the allowable bearing pressure is less than the safe bearing pressure, if 
the footing width exceeds 3 to 4ft. depending on the footing depth. For raft 
and pier foundations it has been suggested(20) that twice the allowable bear- 
ing pressure of footings can be used. For both types of foundations settle- 
ment governs the bearing pressure in all piactical cases and the allowable 
pressure is then given by twice the values estimated from Eq. 7. Since sub- 
mergence of cohesionless soils increases the settlement by about the same 
amount as the bearing capacity is reduced, the allowable bearing pressure 
should be reduced with position of the water table as indicated above for the 
ultimate bearing capacity. 


Bearing Capacity and Settlement of Piled Foundations 


The ultimate bearing capacity of a pile is the sum of the point resistance 

and skin friction and is given by 
Q = (tons) 

where Ap = sectional area (sq. ft.) of pile toe (point), 

A, = surface area (sq. ft.) of pile shaft, 

fs = average unit skin or shaft friction (tons per sq. ft.) 


unit point or toe resistance (tons per sq. ft.). 


and Ip 


In accordance with Eq. 2 the unit point resistance in homogeneous cohesion- 
less soils may be expressed by 
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G = yDN, 
where D = depth of pile point, 


and pile bearing capacity factor depending mainly on the angle of 


internal friction of the soil. 


It has been found(5) that it is safe to take the factor N,, as constant for a 
penetration ratio of D/B 310. For smaller penetrations, the corresponding 
factor can be obtained approximately by linear interpolation between the 
footing factor Ng and the pile factor Np so that 


Gp = yD/Nq + (Np - Ng)D/10B/ for D/B < 10. cca 3 
Similarly it has been indicated that the average skin friction may be ex- 
pressed by 
D 
fs = 6 (11) 
. where § = angle of skin friction, 
I and Kgs = average earth pressure coefficient on the shaft. 


Analysis of field observations on driven piles shows that the shaft friction 
ae K,tan § varies from about 0.25 in loose sand to about 1.0 in dense 
sand. 

From the relationship between N and @ given in Table 1 it can then be 
shown from Eqs. 9 and 11 that approximately 


ND/2 (tons per sq. ft.) 
ND/1000 (tons per sq. ft.). 


Although Eqs. 12 and 13 are in agreement with the experience that both 
unit point resistance and unit skin friction increase with depth and the total 
bearing capacity increases approximately with the square of the depth, field 
conditions are generally too complex for a quantitative estimate on this basis. 
On the other hand, the results of numerous pile loading tests have shown that 
the ultimate bearing capacity of piles in cohesionless soils can be estimated 
from the results of static cone penetration tests with an accuracy sufficient 
for most practical purposes. Thus, it has been found(4, 21) that the observed 
unit point resistance q, of piles varied from about two-thirds to one and one- 
half times the static cone resistance q, and on the average 


dp 
and fs 


Substituting the correlation between static cone and standard penetration 
resistances given by Eq. 1 into 14a, the unit point resistance may be taken as 


4N (tons per sq. ft.) 


average standard (or equivalent) penetration resistance 
(blows per ft.) near pile point. 


dp 
where N 


Similarly field loading tests on driven displacement piles have shown(4) that 
the observed unit skin friction f, of piles varied from about one and one- 
quarter to three times the static skin friction f, on the shaft of a 
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penetrometer and on the average 
f, = 2f... © (15) 


The greater skin friction of piles compared with that observed on a 
penetrometer may be explained by the greater lateral compression of the 
soil during installation of full sized piles compared with that of a penetrome- 
ter, 

To obtain some information about the variation of the skin friction of piles 
with relative density of cohesionless soils, the published skin friction mea- 
surements on driven displacement piles (22,23,24,25) have been analyzed and 
compared with the average static cone resistance measured throughout the 
same depth in a penetration test (Fig. 4). While there is a considerable scat- 
ter of the test results since the skin friction depends not only on the relative 
density of the soil and method of installing the pile, but also on the compressi- 
bility of the soil, pile dimensions and various other factors, it is found that 
on the average 


fs = Gc/200 (tons per sq. ft.) GP OR 
or, substituting Eq. 1, 


fs = N/50 (tons per sq. ft.) ‘chine 


where de = average cone resistance (tons per sq. ft.) 
within depth penetrated by pile, 


and N = average standard (or equivalent) penetration 
resistance (blows per ft.) within same depth. 


Fig. 4 alsc gives the results of a similar analysis of the skin friction ob- 
served in static cone penetration tests.(26,27) The great scatter of the results 
indicates the difficulty of a reliable determination of the skin friction from 
small scale measurements, but on the average these penetrometer friction 
values are about one-half of the skin friction of piles for the reason men- 
tioned above. Hence, the average unit shaft friction on penetrometer is, very 
approximately, 


fe = q-/400 (tons per sq. ft.) 
or, substituting Eq. 1, 
f. = N/100 (tons per sq. ft.). 


On the basis of these correlations the ultimate bearing capacity of driven 
displacement piles can be obtained by substituting Eqs. 14 to 16 into 8, i.e. 


Q = + 2fcAg (tons) 
using the static cone penetration test, 
or Q = 4NAp + NA,/50 (tons) Aen; 


using the standard penetration test results, with an upper limit of the unit 
skin friction of the piles of about 1 ton per sq. ft. For tapered piles Ay may 
be taken as the sectional area at the lower third-point of the embedded pile 
length. 
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Similarly, the ultimate bearing capacity of driven piles with small soil dis- 
placement (e.g. steel H piles) may be taken as similar to that of a penetrome- 
ter with an average unit skin friction 


= fo. 
Hence, after substituting Eqs. 14, 17 and 19 into 8 


Q = GcAp (tons) 
using the static cone penetration test, 


or Q@ = 4NA, + NA,/100 (tons) | 


using the standard penetration test results, with an upper limit of the unit skin 
friction of about 1/2 ton per sq. ft., 


where Ap = net sectional area of toe (e.g. areas of flanges and web), 


and A. = gross surface area of shaft (e.g. all surfaces of flanges 
and web). 


As a check of the proposed relationships between the ultimate bearing 
capacity of piles and the penetration resistance of cohesionless soils, the re- 
sults of pile loading tests(28,29,30) are given in Table 3 together with the 
corresponding estimates using Eqs. 18 and 20. It is found that the observed 
bearing capacities are in fair agreement with the estimated values with a 
maximum difference of about 20 per cent except for the steel H pile. In the 
latter case the estimate is conservative and indicates that the effective toe 
area is somewhat greater than the net section used for the estimate; but fur- 
ther tests are required before an increased toe area could safely be recom- 
mended. 

It should be noted that the correlations between the penetration resistance 
and the point resistance and skin friction of piles have been obtained for satu- 
rated sand. In the infrequent cases where piles rest in dry or moist cohesion- 
less soils the suggested relationships are therefore conservative. On the 
other hand, where the penetration ratio D/B<10, the point resistance has to be 
reduced as indicated by Eq. 10 and is, approximately, 


= G-D/10B (tons per sq. ft.) 
or 4ND/10B (tons per sq. ft.). 


For very small ratios of D/B, Eq. 5 for footings should be used. While the 
ultimate bearing capacity of cast-in-place piles with permanent shell is the 
same as that of displacement piles, a removable shell may reduce the skin 
friction to the lower limit given for steel H piles (Eq. 19), the amount of re- 
duction depending on the degree of permanent displacement of the soil. The 
above equations apply also to piles penetrating through cohesive soil into a 
cohesionless stratum provided the depth of penetration D refers to the co- 
hesionless stratum and the skin friction in the cohesive soil is added to the 
bearing capacity. 

With a factor of safety of 3, the safe pile load is 


Q, = Q/3 
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where Q; is given by Eqs. 18 and 20. The correlations (Figs. 2 and 3) show 
that this factor of safety provides a minimum factor of about 1.8 with respect 
to the most unfavourable results of observed point resistance and observed 
skin friction of displacement piles and penetrometers. 

Since a great part of the bearing capacity of piles in cohesionless soils is 
due to the point resistance, and since piles are generally spaced at not less 
than 2-1/2 to 3 times the pile diameter on centres, the ultimate bearing ca- 
pacity Qe of a pile group may be taken as the sum of the bearing capacity 
Q; of the individual piles, i.e. 


= nQ 


where n = number of piles in group. 


Hence with a factor of safety of 3, the safe load of the pile group is 


Qgs = nQ,/3. 


For piles passing through compressible material into cohesionless soils it 
has been suggested(31) that the safe load should not exceed two-thirds of the 
point resistance ignoring skin friction entirely. 

Within the upper limit of the safe load given by Eq. 24, the allowable load 
on a piled foundation in homogeneous cohesionless soils can be based on simi- 
lar considerations of settlement as the allowable pressure on deep spread 
foundations since the relative density of the soil at a short distance below the 
level of the pile points is sensibly unaffected by the installation of the piles. 

If it is assumed that the load on a pile group is applied to the soil by an area 
located at the elevation of the lower third-point of the length of the piles, the 
allowable load on a small pile group may be taken as that of a pier with a base 
at the above-mentioned elevation and for a large group as that of a raft at that 
elevation. It has also been suggested(20) that in both these cases a pressure 
of twice the corresponding allowable value for a shallow footing of the same 
area may be used and is unlikely to lead to differential settlements exceeding 
3/4 in. On the other hand, the seat of settlement of a piled foundation will al- 
most invariably be in fully saturated soil so that one-half of the loads for dry 
or moist soils have to be used. 

Hence from Eq. 7 since B> 4 ft., the allowable load on a pile may be ex- 
pressed by 


or 


Nd? (1 + 1/B)? (tons) 
12 


for a pile layout in a square pattern and 0.8 Q, for a triangular pattern, 


where d = average pile spacing (ft.) in group, 
B- = width (ft.) of pile group (footing or raft), 


Gc and N = average penetration resistance between depths of 2D/3 and 
2D/3 + B, and other symbols as before. 


It is of interest to note that for an average pile spacing exceeding 3 pile 
diameters, a comparison of Eqs. 22 and 25 indicates that a safe load per pile 
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using a factor of safety of 3 will frequently govern the allowable load and en- 
sure that the maximum settlement of a piled foundation does not exceed about 
2 in. even for a large pile group. This result has been confirmed from an 
analysis of settlement observations of piled foundations in comparison with 
the results of pile loading tests on sands. (32 


CONC LUSIONS 


1) In view of the difficulty and expense of procuring undisturbed samples 
of cohesionless soils, it is frequently most convenient to estimate the relative 
density and bearing capacity of such soils from the results of penetration 
tests on the site. The most important types of tests are standard penetration 
tests, which had previously (1948) been correlated with the relative density 
of sands and the bearing capacity of spread foundations, and static cone pene- 
tration tests, which have for some time (since about 1935) been used for esti- 
mates of the bearing capacity of piles. 

2) Comparison of the results of standard and static penetration tests at 
various locations shows that an approximate relationship exists between the 
standard and static penetration resistances. On the basis of this empirical 
relationship the relative density and approximate angle of internal friction of 
cohesionless soils can be estimated from the static cone resistance. 

3) The suggested relationship between relative density and angle of inter- 
nal friction enables the bearing capacity of spread foundations to be estimated 
from the results of either standard or static penetration tests. Comparison 
of these estimates with the results of plate loading tests on different sites 
indicates that the estimates are somewhat conservative, especially for small 
footing widths. 

4) From the empirical correlation between standard and static penetration 
tests and an analysis of the skin friction measured in pile loading tests and 
static penetrometers, a method is developed to estimate the bearing capacity 
of driven displacement piles and piles with small soil displacement (for in- 
stance, steel H piies) from the standard penetration resistance. Comparison 
of these estimates with the results of some pile loading tests shows fair 
agreement, except for the steel H pile, where the effective toe area is some- 
what greater than the net section used for the estimate. 

5) The proposed procedure can also be used to estimate the allowable load 
on spread and piled foundations in cohesionless soils from the results of 
standard and static penetration tests. 
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Details of Penetration Tests 
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FIG.3. COMPARISON OF STATIC CONE RESISTANCE AND RESULTS 
OF STANDARD PENETRATION TESTS 
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TABLE 1 
RELATIONSHIP BETWEEN RELATIVE DENSITY, PENETRATION RESISTANCE 


AND ANGLE OF INTERNAL FRICTION OF COHESIONLESS SOILS 
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TRIAXIAL SHEAR TESTS ON PERVIOUS GRAVELLY SOILS 


Wesley G. Holtz,! M. ASCE and Harold J. Gibbs,? A.M. ASCE 
(Proc. Paper 867) 


SYNOPSIS 


This paper was prepared to present the results of recent tests to deter- 
mine the shear characteristics of free-craining gravelly soils. These tests 
are part of a research program for determining the shear characteristics of 
all types of gravelly soils which is in turn only a portion of a larger research 
program designed to study all of the physical properties of soils containing 
gravel particles. Large-size triaxial shear test equipment was fabricated 
for making the shear tests. The testing program was formulated to provide 
information on the relations between shear resistance and (1) density, (2) 
amount of gravel, (3) gradation, (4) maximum particle size, and (5) particle 
shape. Also included in the paper is a discussion of the effect of testing 
speed on shear test data and the effect of specimen size for any maximum 
particle size up to 3 inches. The conclusions reported in the paper are based 
on an analysis of 183 individual shear specimens. 


INTRODUCTION 


In the past, laboratory soil tests have been largely limited to the finer 
fraction in which gravel particles have been removed. As a general rule, 
tests are seldom made on soils with particles larger than 3/16 inch (No. 4) 
size. This has been largely due to the fact that larger specimens are re- 
quired if the coarse fractions are used and the testing costs are increased 
greatly. The Bureau of Reclamation has recognized a need for information 
on the physical properties of the total fractions of materials containing 
coarse particles because so many of the materials being used in our earth 
dams, earth canal embankments, and earth canal lining works contain these 
coarse particles. It is often good engineering practice to use gravelly soils 


Note: Discussicn open until May 1, 1956. Paper 867 is part of the copyrighted Journal 
of the Soil Mechanics and Foundations Division of the American Society of Civil En- 
gineers, Vol. 82 No. SM 1, January, 1956. 

1. Chief, Earth Laboratory Branch, U. S. Bureau of Reclamation, Denver, 
Colo. 

2. Head, Special Investigation Section, Earth Laboratory Branch, U. S. 
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when practicable because of the high stability characteristics of these materi- 
als. 

For these reasons a broad research program was outlined by which all the 
physical properties of the various types of gravelly soils would be studied. 
This broad program is being pursued as time permits. Some research on 
compaction, (1) permeability, and consolidation, (2,3) in addition to the shear 
results presented herein have been completed. These shear test results were 
obtained from a rather extensive study on free-draining sand-gravel mixtures. 
Tests have now been started to study the shear characteristics of gravelly 
soils of low permeability. 

In the series of tests on free-draining sand-gravel mixtures reported here- 
in, an attempt was made to determine the relations between shear resistance 
and (1) density, (2) amount of gravel, (3) gradation, (4) maximum particle size 
and (5) particle shape. In order to properly study these factors, it was also 
necessary to investigate the effect of speed of testing on shear test results for 
the various specimen sizes used and the effect of specimen size on shear test 
results for the various maximum particle sizes used. This information is also 
covered in this paper. In this entire research study 31 complete shear tests 
were made. This involved the testing of 183 individual shear specimens which 
is an average of about six specimens per shear test. Figure 1 is a line chart 
which shows the complete testing program. 

For the sake of brevity much of the detailed test information has not been 
included, but the over-all results are discussed and appropriate examples 
given. The principal details of the testing program are contained in the 
Bureau of Reclamation publication EM-350, “First Progress Report of Re- 
search on Triaxial Shear Testing of Gravelly Soils,” dated July 17, 1953. The 
placement density was controlled on the basis of relative density which has an 
advantage from a practical standpoint. Therefore, the data are also presented 
from an “initial placement” relative density standpoint. Actually, in the an- 
alysis of the data, load-void ratio curves were plotted from the test data and 
the results were also studied from the standpoint of failure conditions. It was 
decided that the “initial placement” conditions provided a suitable level for 
the comparisons presented herein. 

The plots of the Mohr envelopes of limiting shear resistance were not ex- 
actly straight lines but were generally slightly curved convex upward. From 
a research standpoint straight line determinations would be inadequate and 
sometimes misleading because slight variations in curvature would cause 
fluctuations in computed straight line values; therefore, plots of the envelopes 
are shown for comparing results. Straight line interpretations are often used 
in design work and are also used herein for discussion purposes because of 
their simplification. When so used herein they are called computed values, 
being computed from the most reasonable straight line envelope to the shear 
circles. The pertinent test data for all specimens tested have been recorded 
in Table I. 


Equipment 


The Bureau of Reclamation Earth Laboratory has (1) small triaxial shear 
equipment for testing 1.375-inch diameter by 3-inch long specimens, (2) medi- 
um size equipment for testing 2.625-inch by 7-inch and 3.25-inch by 8.125 or 
9-inch specimens, and (3) large equipment for testing 6-inch by 15-inch and 
9-inch by 22.5-inch specimens. The details of these shear machines have 
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been described elsewhere, (4,5) and the large machine will be discussed here 
only in a brief manner as needed to understand the data presented. The three 
machines operate in the same general manner and permit determinations of 
the same over-all data. The two smaller machines contain axial loading de- 
vices built as an integral part of the machine, whereas the large equipment 
utilizes a compression testing machine. 

Figure 2 shows the large shear testing equipment. The compression ma- 
chine for applying the axial load is equipped with automatic stress-strain re- 
corder and rate indicators. The pressure chamber was designed for a maxi- 
mum working lateral pressure of 100 psi. It consists of aluminum top and 
bottom plates, aluminum or clear plastic chamber cylinders of 14-inch I.D., 
and six 1-inch stainless steel bolts to hold the assembly together. The piston 
through which the axiai load is applied to the specimen was made of stainless 
steel 2 inches in diameter. The piston passes through a bushing in the top 
chamber plate, and “O” rings recessed into the bushing provide a positive 
seal around the piston. A small stainless steel plate attached to the top of 
the piston provides a bearing plate for the head of the compression machine. 
A dial indicator of 1-inch travel is attached to this plate. The indicator 
travels along a slanting bar attached to the top chamber as the specimen is 
compressed. By this means axial deformation can be measured for as much 
as 7 inches or about 30 percent axial strain. 

The specimen end plates were made from 1-inch aluminum plate and re- 
cessed on the specimen side for inserting porous bronze or perforated plates 
0.125 inch thick. The top and bottom end plates are fitted with Saran tubing 
and brass cocks for connecting to drainage or saturation tanks and/or pore 
pressure cells. Saran tubing was used because of its flexibility and ability to 
stand 100 psi outside pressure without collapsing. The rubber membranes 
for enclosing the soil specimen are 6 inches in diameter, 16.5 inches long, 
and 0.03 to 0.06 inch thick for the 6-inch diameter specimens and 9 inches in 
diameter, 24 inches long, and from 0.03 to 0.06 inch thick for the 9-inch di- 
ameter specimens. Stainless steel “Aero-Seal” clamps of the worm screw 
type are used to clamp the membrane to the end plates. 

The lateral pressure is applied by means of air pressure on the water- 
filled chamber. The lateral pressure, and pore pressures at each end of the 
specimen, are measured either by no-flow pressure cells or automatically 
through pressure transducers by means of an electronic strip recorder. 
Water gages with Vernier reading glasses are mounted on the back side of the 
control panel for reading the volume change of the specimen by registering 
the amount of water expelled or taken up by the chamber during the test. 

The equipment allows the performance of Unconsolidated-Undrained, 
Consolidated-Undrained, or Drained tests. All of the tests on the free- 
draining material discussed herein were of the latter type, and were saturat- 
ed prior to shearing. Drainage was accomplished through the bottom end 
plate and pore pressure measurements were taken from the top end plate to 
insure that complete drainage (zero pore pressure) conditions existed as each 
increment of load was applied. The rate of strain necessary to produce this 
condition was not exceeded. The above facilities allow the measurement of 
the following data during the test: axial load, axial strain, volume of speci- 
men, volume change of specimen, permeability, volume of drained water, and 
pore pressure. 
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Materials 


The gradation of the various soil and gravel materials initially selected 
for study are shown in Figure 3a. The soil was a pervious, medium-to- 
coarse, sand consisting predominately of feldspar and quartz particles. To 
this was added varying amounts of 3/16-inch to 3/4-inch, and 1-1/2-inch, 

and 3-inch gravel to obtain the gradation of the curves shown, The gravel and 
coarse sand particles were predominately gneiss, granite and schist. The 
sand and gravel was obtained from a river deposit and the particles were of 
subangular to subrounded shape. The specific gravity of the sand and sand- 
gravel mixtures varied from 2.52 to 2.63. As the material was reused numer- 
ous times, a gradation test was made after completion of the tests to check 
particle breakdown. It was found that particle size changes were very minor. 
The second group of materials were tested in connection with the design 

of a rockfill structure. These materials consisted of very sharp angular 
particles of porphyritic quartz monzonite and diabase obtained from rock 
quarry operations. The gradations of the materials tested are shown by the 
curves in Figure 3b. The particles of this material tended to break down ap- 
preciably upon reuse. The specific gravity of the finer and coarser grada- 
tions was 2.68 and 2.85, respectively. 


Specimen Preparation 


The shear specimens were prepared by compacting moist material to the 
desired density in a hinged three-section mold. Care was taken to insure 
that segregation of the material did not occur. After the specimen had dried 
slightly in the mold, it was removed and the membrane expanded and slipped 
over the specimen. In some instances where the gravel-to-sand proportion 
was high, the specimen was compacted in the membrane and a vacuum applied 
to the end plates so that sufficient strength was obtained for handling. 

Control of placement was based upon relative density for each material or 
mixture. The densest state (100 percent relative density) was determined by 
standard laboratory procedure(4) which consisted of compacting saturated 
material into a 0.5-cubic-foot container by intense vibration. The loosest 
state (0 percent relative density) was determined by standard laboratory pro- 
cedure(4) which consists of carefully hand placing the material in a 0.5-cubic- 
foot container. The desired specimen density was then computed on the basis 
of the following equation: 


y= min 


Ymin 
1-D, (1 - Ymax) 


100 


where 

Y = the density in densest state 
» AOS... the density in loosest state 

Y =the desired specimen density 

D. = the desired relative density in percent 


Specimens of the river sands and gravels were placed at 50 and 70 percent 
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relative density to obtain data for observing the effect of density on shear 
strength. These ranges of density were considered to be within the range of 
economical and practical placement conditions. The specimens of the quarry 
rock were placed at 50 to 90 percent relative to observe the effect of extreme 
density ranges. 

Placement moisture was not a controlling factor in this study since the 
materials were free draining and non plastic. The specimens were prepared 
with the material in a moist condition and were saturated by percolating water 
through them after application of the chamber pressure but prior to loading 
axially. 


Studies on Rate of Load Application 


Prior to conducting the planned group of tests on the large shear specimens, 
it was necessary to know what effect the specimen size had upon the shear 
values obtained. In order to study this factor properly it was also necessary 
to perform some tests which would provide information on the effect of testing 
speeds. The study on rates of loading included three different rates which 
were classified as “fast,” “medium,” and “slow.” These represented rates 

of axial strain of approximately 1.8, 0.36, and 0.086 percent per minute, re- 
spectively. Figure 1, previously referred to, shows the group of shear tests 
which were made to determine this factor. Fourteen complete shear tests of 
the river sand and gravel were analyzed for this phase of the work. The dif- 
ferent rates of testing were tried for six conditions in which the density, 
specimen size, and sand-gravel content were varied. 

It became evident early in the program that the rate of testing had practi- 
cally no effect upon the shear test results for all of the specimen sizes. In 
other words, these materials were sufficiently free draining to allow all the 
stresses and strains to take place rapidly within the range of testing periods 
tried. Figure 4, a and b, shows the shear envelopes obtained in this phase of 
the study. In these figures, curves for similar material, and placement den- 
sity and specimen size conditions practically coincided. There was no con.. 
sistent trend in the slight differences which did exist. Therefore, the medium 
rate was used for the remaining study because it was most practical from a 
test operation standpoint. 

Volume change information taken during these tests were plotted in the 
form of void-ratio versus deviator stress curves for the purpose of check of 
testing trends in addition to those shown by the shear envelopes. These 
curves also showed practically no effect from variations in rate of loading. 


Studies on Specimen Size 


Four sizes of specimens were selected for studying the effect of specimen 
size on measured shear strength. The diameter and length dimensions were 
1.375 by 3 inches, 3.25 by 8.125 inches, 6 by 15 inches, and 9 by 22.5 inches 
and for purposes of discussion are designated as very small, small, medium 
and large, respectively. These specimens all had length/diameter ratios of 
2.5 except the very small size which had a length/diameter ratio of 2.18. Six- 
teen complete shear tests on the river sand and gravel were analyzed for this 
phase of work. Different size specimens were tested for six conditions in 
which density, rate of testing, and sand-gravel content were varied. 

The results of the “specimen size” tests are summarized on Figure 5, a, 
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b,c. Part “a” is significant because it shows a comparison of test results for 
all specimen sizes when the material tested consisted only of the sand frac- 
tion (passing No. 4 seive size) placed at 70 percent relative density. The very 
Small specimen size was the only one showing marked difference in strength. 
It is felt that this size is too small for material having appreciable coarse 
sand of 3/16-inch maximum size. Part “b” of Figure 5 shows the results of 
tests on the small, medium, and large specimens with material of 80 percent 
sand and 20 percent gravel of 3/4-inch maximum size. The results of these 
tests were practically identical. Part “c” of Figure 5 shows the results of 
similar tests where the gravel content was increased to 50 percent. Here it 
can be seen that gravel of this amount significantly affected the results and 

a higher shearing resistance was indicated for the small size specimen than 
the medium and large sizes which gave near similar results. The apparent 
conclusion is that the small 3.25-inch diameter) specimen is not large enough 
for 3/4-inch maximum size material when appreciable amounts (over 20 per- 
cent) of coarse particles are present. For this reason, only the 6-inch and 
9-inch diameter specimens were used for the remainder of the study on sand- 
gravel mixtures. 


Effect of Relative Density 


Figures 4, a, b, and c are plots of the shear envelopes for tests on large 
specimens which show the effect of relative density on the shear strength. 
Twelve complete shear tests on large specimens were analyzed for this phase 
of study. Different relative densities were tried for four conditions in which 
gravel content, particle size, and particle shape were varied. 

In all instances, the envelopes for the higher density specimens were ap- 
preciably higher than those of the lower density specimens. The computed 
friction values (tan@) for the river sand placed at 70 percent and 50 percent 
relative density were 0.75 and 0.65, respectively. The computed friction 
values for the river deposit material consisting of 80 percent sand and 20 per- 
cent gravel to 3/4-inch size were 0.78 when placed at 70 percent relative den- 
sity and 0.69 when placed at 50 percent relative density. The finer quarry 
material, which consisted of 35 percent sand size and fines, and 65 percent 
gravel to 3-inch size, showed computed friction values of 0.85 when placed at 
70 percent relative density and 0.80 when placed at 50 percent relative density. 
The coarser quarry material which consisted of about 18 percent sand size 
and fines, and 82 percent gravel to 3-inch size showed computed friction 
values of 0.98 at 90 percent relative density and 0.85 at 65 percent relative 
density. 


Effect of Gravel Content 


The effect of gravel content on the shear strength of earth material was 
the principal factor of this study although the other factors studied are, of 
course, closely related. This factor was studied by analyzing thirteen com- 
plete shear tests. The typical variation in shear strength, with gravel con- 
tents varying from zero to 65 percent, are shown by the envelope plots, Fig- 
ure 6. Part “a” of the figure shows envelopes for zero, 20 percent, 35 per- 
cent, 50 percent, and 65 percent gravel for the river sands and gravels of 
3/4-inch maximum size tested at 70 percent relative density by large speci- 
mens. The computed friction values (tan @) were 0.72, 0.76, 0.81, 0.84, and 
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0.78, respectively. Other test results for (1) zero, 20 percent, and 50 percent 
of 3/4-inch maximum size river gravel contents at 70 percent relative density 
and tested by medium-size specimens, (2) 20 percent and 50 percent of 1-1/2- 
inch marimum-size river gravel, and (3) 20 percent and 50 percent of 3-inch 
maximum size river gravel, showed similarly increased friction values with 
increased gravel contents, up to 50 percent. 

Part “b” of Figure 6 shows the envelopes for two of the quarry materials 
which had 65 percent and 80 percent grave: contents, with a maximum particle 
size of 3 inches. For the series of specimens placed at 50 percent relative 
density for the finer and coarser materials, the computed friction values were 
0.80 and 0.85, respectively. 

It can be generally stated that the friction values increase, other factors 
being equal, as the gravel contents increase up to gravel contents of, say, 50 
to 60 percent depending on the maximum particles sizes. When gravel con- 
tents of this amount are exceeded, little or no increase in shear strength is 
evident; in fact, some reduction in strength may occur. This is probably 
caused by the material becoming less well-graded. It will also be noted 
(Table I) that the actual densities, for the 3/4-inch maximum size river grav- 
el mixtures placed at 70 percent relative density, increased up to about 50 
percent gravel content, while the actual density did not increase for the 65- 
percent gravel content. 


Effect of Maximum Particle Size 


The other principal factor studied in these experiments was the effect of 
maximum particle size upon shear strength. The results of six complete 
shear tests on large specimens were analyzed for this phase of study. Fig- 
ure 7 is a plot of shear envelopes to demonstrate this factor. Part “a” of the 
figure shows envelope plots for shear tests on large specimens with 20 per- 
cent river gravel content and placed at 70 percent relative density. The en- 
velopes for 3/4-, 1-1/2-, and 3-inch maximum particle size are practically 
identical although the straight line computed friction values (tan@) varied 
slightly, being 0.76, 0.78, and 0.79, respectively. Actually, in the case of 
these 20-percent gravel contents the larger gravel pieces represented only a 
very small portion of the total material in number of particles. Part “b” of 
the figure demonstrates the envelopes for 3/4-, 1-1/2-, and 3-inch maximum 
particle sizes with river gravel contents of 50 percent and placed at 70 per- 
cent relative density. Here again the envelopes are almost identical although 
the computed shear values were 0.84, 0.84, and 0.88, respectively. No com- 
parisons of maximum particle size was made on the quarry materials. From 
these studies, it is our opinion that the maximum gravel particle size from 
3/4- to 3-inch has little effect upon shear resistances because the larger 
heavy particles are actually few in number, other conditions being equal. 
There may be a slight increase in strength with increase in maximum size, 
particularly for gradations where the larger pieces result in a greater num- 
ber of particles. 


Effect of Particle Shape 
Although the program as initially set up with the river gravels did not en- 


compass a study of particle shape effect on shear strength, a recent design 
study with the quarry materials provided data for a few brief comparisons. 
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One series of the river sand and gravel materials were similar to one series 
of the finer quarry materials in that both series consisted of specimens of 65 
percent of gravel to 3-inch maximum particle size and placed at 70 percent 
relative density. As shown by the shear envelopes of Figure 8, the shear 
strength of the more angular material is somewhat higher than the shear 
strength of the subangular to subrounded material. The computed friction 
values (tang) for the materials, under the conditions described above, were 
0.78 and 0.85. It was found, moreover, that as the density of the coarser 
sharp quarry material was increased to 90 percent relative density, very 
high friction values (tan ¢ of 0.98) could be obtained. 


CONCLUSIONS 


The conclusions listed below may be stated as a result of the experiments 
reported herein on the river gravel deposit and quarry materials. Figure 9 
is a summary plot of the test results obtained and shows the range of materi- 
als studied and the range of shear values obtained. 


a) For the free-draining sand and sand-gravel mixtures used in this study, 
variation of the rates of axial-strain from 1.8 to 0.086 percent per minute 
had no significant effect upon the shear strength. An intermediate rate was 
selected as the practical from a test operation standpoint. 

b) Four sizes of shear specimens were studied varying from 1.375 to 9 
inches in diameter. It was concluded that the 1.375-inch diameter size was 
too small for tests on the minus 3/16-inch fraction if appreciable amounts of 
coarse sand were present and that the 3.25-inch diameter specimen was too 
small for 3/4-inch maximum size material when appreciable (over 20 percent) 
gravel sizes were present. It is believed that 6-inch diameter specimens are 
suitable for materials having a maximum particle size up to 1-1/2 inches, 
and the 9-inch diameter is suitable for materials, of normal gradation, having 
up to 3-inch maximum size particles. 

c) The “initial placement” relative density of the sand or sand-gravel mix- 
tures had appreciable effect on the shear strength. The lowest computed 
friction value (tan ¢ = 0.65) was obtained on the river sand material placed 
at 50 percent relative density while the highest value (tan ¢ = 0.98) was ob- 
tained on the quarry material having 82 percent gravel sizes and placed at 
90 percent relative density. 

d) While the size of the gravel particles above 3/4 inch itself had no ap- 
preciable effect on shear strength because they were few in number, the 
amount of gravel in the mixture had a very significant effect. The shear 
strength appears to increase appreciably as the gravel content is increased 
up to 50 to 60 percent, depending on the maximum size. After this point is 
reached the material becomes less well-graded and the actual density does 
not increase (for constant relative density) and likewise the shearing resist- 
ance does not increase, or even may become less. 

e) The shape of the particles affected the friction of the material signifi- 
cantly. An increase in shear strength was observed for the angular quarry 
material over the subrounded to subangular river materials and density ap- 
peared to have a particularly significant effect on the friction of the angular 
materials. 
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FIGURE 2 is 
TRIAXIAL SHEAR EQUIPMENT FOR LARGE SIZE SPECIMENS . 
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Test B--508 relative density-- 
fast rate 

Test D--50% relative density-- 
slow rate 

Test J--70% relative density-- 
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SHEAR STRESS--PSI 
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FIGURE 6 - EFFECT OF GRAVEL CONTENT 
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FIGURE 7 - EFFECT OF MAXIMUM PARTICLE SIZE 
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Table 1 
| SUMMARY OF TEST DATA Sheet 1 of 4 |. | 
~ § > enve. 
| | beh | 2. | computation HA 
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A 3/16 36 103.2 3.3 13.4 | 108.5 |0.570 
39 103.8 6.4 21.5 | 105.7 | .552 
ho | 105.3 12.7 48.1 | 109.2 | .503 
104.2 25.0 107.6 527 
ha | 108.5 50.0 | 161.7 | 112.2 | 
43 104.9 9.9 308.6 | 113.4 | 
Ay (2 
3/16 | 13.6 | 105.5 | .555 
915 | 10k.2 20.6 | 106.8 | 
916 | 104.8 32.9 | 105.9 | .5k8 
917 | 106.8 60.9 | 109.2 | .502 
29 | 105. | 
Avg ims he 
le o| 3/16 be 104.7 3.1 1.9 | 105.4 | .555 
45 104.6 6.2 20.4 | 105.6 | 540 
103.1 12.6 3.2 | 109.1 | .537 
| 25.2 &.8 | 107.2 | .530 
| 108.2 49.9 | 166.2 | | 
4 4g 99.9 301.5 | 14.9 | 
Avg 104.0 
D 3/16 ‘Siow 920 3.5 | 12.5 | 105.3 | .560 = 
6.0 18.8 | 105.6 
gee 12.3 36.3 | 105.1 | 
923 2.9 &.9 | 108.0 | .520 
gek 49.9 127.8 108.3 -516 = 
925 99-9 | 259.2 | 1.6 | 
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17 50.1 | 20 
a | 100.0 35 6 | 
3/16 31 4.3 | .500 
| 7.2 +506 
100.2 | 29% +469 
Ave | 
3/6 605 | 3.5 +500 
609 | 6.7 x -503 : 
60 | 12.9 b | +539 
25.8 9 +460 
607 50.3 156 | 
604 100.3 29) 
Ave -T2 &.7 
3/6 905 3-9 ay -522 
906 6.7 2 -509 
901 25.4 91] 
goe 50.1 | 16d 
903 100.2 | 308 
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3.3 23.7 | | 
6.3 | 32.8 | 10.6 | 
12.6 | 65.6 | 112.2 | 
25.0 | 120.6 | 112.2 | 
69.9 | 196.7 | 106.6 | 
99.9 | 338.2 | 125.2 | 
10.489 0.80 4.7 
o| 3/16 Large 907 110.6 17.1 
913 108.2 25.9 
909 108.7 52.8 
ave 109.5 3.7 
20 | 926 111.6 3.2 20.6 | 122.6 | .¥73 : 
927 110.6 6.3 35.5 | 211.7 | .493 
928 211.5 12.6 | 46.9 | 121.5 | 
ge9 111.9 25.2 66.4 12.4 455 
930 no.4 $0.2 | 156.2 | 113.2 | 
112.1 100.0 | 273.9 | 116.2 | 
933 109.9 an 
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20 | 939 3.6 | 2.0 | 125.3 
2 12.8 | 63.5 | 114.8 
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a ~ > enve lope 
610 117.6 | 0.392 
6 117.1 .397 
6.7 17.1 | .397 
613 115.3 | 
15 126. lg. 373 
3/6 33 us.6 3.3 2.6 | 125.0 | 
3 115.3 6.6 | 125.1 | 
ng 5.5 12.8 63.1 | 116.2 | 
320 1.8 25.2 %.5 | 116.0 | 
322 5.6 50.3 177.6 | 117.6 | .392 
322 115. 100.1 328.2 | 19.4 | .372 
959 19.0 3.5 26.2 | 18.9 | .376 
960 120.2 6.5 | 119.7 | .360 
961 119.0 12.9 73-5 | 119.9 | .358 an 
962 119.1 25.4 120.6 | 120.2 | .353 
963 19.3 205. | 120.7 | 
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2 ve 3 
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H 623 100.4 3n.5 285 
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308 125.8 50.3 261.1 | 126.1 | .267 
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*These test specimens contained such a large proportion of coarse angular particles that large voids existed on 


the surfaces of the specimens. 


Therefore, it was not possible to obtain accurate volume change data. 
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REVIEW OF THE PAS’ YEAR 


On Friday, October 28, at the Society’s Annual Convention in New York, 

the Division Executive Committee held its regular annual meeting. The 
Annual Report of the Division was presented, discussed, and accepted; and 
several additional matters received attention. A summary of some noteworthy 
features of the Annual Report and the Executive Committee Meeting follows: 

Committee on Earth Dams: T. M. Leps was appointed Chairman to fill the 
vacancy created by the resignation of former Chairman W. G. Holtz. Mr. 
Holtz will carry on as a member of the Committee. The Committee held a 
meeting on October 28 at the Society’s Annual Convention in New York. 

Committee on Grouting (R. E. Davis, Chairman): The Task Committee on 
Chemical Grouting has completed a report containing, in addition to a brief 
statement of principles and procedures, summarized accounts of about ninety 
actual field applications of chemical grouting, a list of patent abstracts dating 
from 1890, and an extensive bibliography. The Executive Committee will con- 
sider formal publication of this report if sufficient interest is expressed by 
the Division membership. Anyone interested in securing a copy in the event 
of publication is therefore requested to communicate with the Division Secre- 
tary. 

The Task Committee on Cement Grouting has prepared a preliminary re- 
port that is at present under discussion by the Committee. A. W. Simonds 
has replaced V. L. Minear as Chairman of this Committee. The report of the 
Task Committee on Clay Grouting is in draft form. 

Committee on Engineering Geology (K. Terzaghi, Chairman): All task com- 
mittees have been active in the preparation of symposia and drafts of their 
work. The major efforts of the Task Committee on Products of Weathering of 
Bedrock and their Engineering Properties were directed toward preparing for 
a symposium conducted at the Annual Meeting of the Geological Society of 
America at New Orleans on November 7, 1955. The Task Committee on the 
Influence of Geological Factors on Tunnel Construction presented a symposium 
at the ASCE Meeting in St. Louis on June 16, 1955. Meanwhile, the Chairman 
has prepared a paper on the influence of geological factors on tunnel construc- 
tion. This paper is in the hands of the Committee at present. 

Committee on Publications: As reported in the September Newsletter, 

H. B. Seed will assume the chairmanship of this Committee, to succeed 
Hudson Matlock, on January 1, 1956. The Committee is engaged in the devel- 
opment of an effective procedure for handling one of the Divisions most ardu- 
ous and demanding functions. The Committee reported partial or complete 
review of 50 papers during the year that ended in October. Twenty-eight of 
these were still in process of review; and of the remainder, six had been 
recommended for publication, one had been referred to another division, eight 
had been returned to their authors for suggested modification, and seven had 
been declined. 

Committee on Cooperation with Local Sections (R. F Klauck, Chairman): 

In view of the light response to the efforts of this Committee, the Fxecutive 
Committee decided to disband it. An expression of the appreciation and thanks 
of the Executive Committee on behalf of the Division will be conveyed to Mr. 
Klauck and his Committee for their services. The Newsletter will attempt to 
carry on the functions of this Committee. 
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Committee on Glossary of Terms and Definitions (R. E. Fadum, Chairman): 
Several meetings have now been held with the object of reconciling differences ae ; 
between the proposed recommendations of this Committee and those of an a 
ASTM committee having parallel objectives. 
Committee on Road and Airfield Soil Problems (C. R. Foster, Chairman): a 
The members of the Committee, as well as of a board of consultants, have + 
been selected; and preparations for beginning work are veing made. aoe 
Joint Committee of the U. S. National Council and the Soil Mechanics 
Division (D. W. Taylor, Chairman): The Committee membership is complete. 
There has been no activity to date. .? 
Committee on Technical Sessions (S. D. Wilson, Chairman): The results et 
of this Committee’s efforts have been impressively in evidence at San Diego, . | 
St. Louis, and New York. 


CURRENT RESEARCH AND DEVELOPMENT A 


With renewed expression of appreciation, not to mention a sinking heart, 
we present the last installment of the material contributed more than a year 
ago by W. J. Turnbull and his associates at the Waterways Experiment Station 
in Vicksburg, Mississippi. After wallowing for approximately fifteen months 
in the luxury of a continuous backlog, we look with pampered horror on the 
prospect of having to get out on the road and hustle again. 


Report from the Waterways Experiment Station 


Soil Classification. A study of soil classification has been conducted over Re. 
several years by the Office, Chief of Engineers for use in the Corps of Engi- a 
neers. This study has resulted in adoption of the Unified Classification Sys- - 
tem based on Casagrande’s Airfield Classification System. A manual has 
been published by the Waterways Experiment Station explaining the use of the a 
unified system of soil classification with appendices describing characteris- po 
tics pertinent to embankments, foundations, roads and airfields. Additional 
studies are in progress which may lead to the refinement of the present 
system. 

Pipe Cover Requirements for Airfields. The Flexible Pavement Labora- 
tory of the Waterways Experiment Station is making a study of minimum cov- 


er requirements for pipe under military airfield pavements, with emphasis Psi | 

on the requirements for high-pressure tires and multiple wheel assemblies. of f 

Requirements for both rigid and flexible pipes are being studied. ys 


Rational Design of Flexible Pavements. A series of related studies is be- 
ing conducted by the Flexible Pavement Laboratory of the Waterways Experi- 
ment Station to develop a more rational method for designing flexible airfield 
pavements than the empirical California Bearing Ratio method currently being 
used. The approach is being made through measurement of soil stresses and 
strains in field test sections to develop means of computing the stresses in- 
duced by load and through laboratory studies to develop means of measuring 
the true ability of soils and base courses to resist the stresses. re 

Test sections have been constructed of homogeneous silty clay and homo- 
geneous sand, and measurements of soil pressures and strains made under 
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conditions of load simulating airplane wheel configurations. Pressures are 
measured by earth pressure cells; strains by means of specially constructed 
deflection gages. Limited measurements of shear stresses have been made 
with a specially constructed shear cell. It is planned to extend the test sec- 
tion work to layered systems. Theoretical computations have been in vrog- 
ress concurrently to evaluate stresses and strains in homogeneous and lay- 
ered systems. Laboratory tests of the triaxial type are being conducted in 
the laboratory in an effort to duplicate the stress-strain systems measured 
in the field test sections. 

Some developmental work is being done on earth pressure cells, deflection 
gages, and shear cells in order to provide the best possible instruments for 
field measurements. A large-diameter triaxial shear device, 36 in. diam. by 
6 ft. high, has been constructed and is being used for the study of the action 
of soil around measuring instruments. 

Soil Compaction Studies. The Flexible Pavement Laboratory of the Water- 
ways Experiment Station has been conducting a study of soil compaction for 
the Office, Chief of Engineers for over ten years. The more recent work in 
this study is described below. Field test sections have been constructed in 
which the soil types, moisture contents, and coverages of compaction equip- 
ment have been carefully controlled. Studies have been made to determine 
the effect of varying the size of feet on a sheeps foot roller from 7 to 21 sq. 
in. and compacting a silty subgrade. Studies are in progress on the effect of 
(1) varying the tire pressure of a heavy rubber-tired roller from 50 to 150 
psi and compacting a silty clay subgrade; (2) varying the thickness of com- 
pacted lifts from 6 to 18 in. and compacting a silty clay subgrade with rubber- 
tired rollers having tire pressures of 90 and 150 psi; (3) compacting a crushed 
limestone base course with a 90-psi rubber-tired roller and varying the water 
content, size of aggregate, and number of coverages of the roller. 

Office studies also are being made to determine the compaction require- 
ments for base course and subgrade materials for flexible pavements for 
military airfields. A general review is being made of available airfield and 
test section data to study general compaction requirements and procedures 
with particular reference to the effect of yielding subgrades on the compaction 
of overlying materials. 

Landing Mat Studies. The purpose of the investigation being made by the 
Waterways Experiment Station is to provide design criteria which will indi- 
cate the service life of airplane landing mats in terms of wheel loads, tire 
pressures, and subgrade requirements. The investigation consists of both 
empirical and theoretical approaches. The empirical approach consists of 
traffic tests conducted on various types of landing mats laid over a range of 
typical subgrade materials, with and without base courses, using a range of 
single and multiple airplane wheel loads. Data from the tests are used to 
prepare California Bearing Ratio design curves for landing mats for various 
categories of military operations. 

The theoretical approach includes studies to determine the physical char- 
acteristics of the mat, analytical studies of stresses in the mat, model tests 
of mat behavior, and stress distribution studies on the combined behavior of 
the mat and the underlying soil. 
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1956-1 Soil Mechanics and Foundations Division 
CONTRIBUTIONS TO THE NEWSLETTER 


It has recently come to our attention that at least one Newsletter reader 
outside the United States has been under the impression that contributions 
from countries other than the United States are by specific invitation. This is 
not the case; and we take this opportunity to clarify any statement of ours that 
may have ied to such a misunderstanding. 

We occasionally find it advisable to do a little pump-priming; but the oper- 
ation is never performed on a grand scale, nor with any more complicated ob- 
jective than to get something started. Any item of news or other information 
that contributes to the development of general activity and interest in soil me- 
chanics and foundation engineering is a valid entry for the Newsletter. This 
is one publication that welcomes unsolicited manuscripts. 


“SOIL MECHANICS, PURE AND SIMPLE” 


Thus did Mr. Mason Lockwood, of Lockwood and Andrews, Houston, Texas, 
title the talk he gave at the Soil Mechanics Luncheon at the Annual Convention 
in New York. He spoke as an engineer of long experience, who, while not a 
specialist in soil mechanics, can remember its birth, has observed its growth 
with interest, and has had frequent occasion to use what it has to offer. Fol- 
lowing a brief outline of the history of modern soil mechanics and a diverting 
account of some of his own early contacts with the sprightly newcomer when 
both of them were wearing their working clothes, he reached what we took to 
be his salient point, namely, the recommendation that specialists in soil me- 
chanics begin making a conscious effort to remove some of the mystery that 
still surrounds it in the eyes of many non-specialists—an effort, that is, to 
make soil mechanics pure and simple. 


ANNOUNCEMENTS 
DALLAS CONVENTION 
Hotel Baker February 13-17, 1956 


John A. Focht, Jr., of the Division’s Technical Sessions Committee, has 
been in charge of the arrangement of the soil mechanics program for this 
meeting. Plans have been completed for two technical sessions, one devoted 
primarily to some aspects of consolidation and settlement, and the other to 
a variety of pile loading tests and related matters. A luncheon has also been 
planned, complete with presiding officer (Raymond Dawson), speaker (Raymond 
Mason), subject (Problems of Earth Dam Design in Rain Forests), slides 
(color), and, we presume, something to eat. Contributions of technical papers 
will include C. I. Mansur, R. I. Kaufman, W. G. Shockley, and T. B. Goode, all 
of the Waterways Experiment Station; J. B. Eustis, Eustis Engineering Com- 
pany, New Orleans; F. S. Mitchell, Jr., of the Texas Highway Department, 
Houston; L. A. Du Bose, A. and M. College of Texas; F. G. Bryant, University 
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of Texas; Bramlette McCleiland, of McClelland Engineers, and Lindsey Lips- 
comb, of Humble Oil and Refining Company, Houston. Presiding at the tech- 
nical sessions will be W. J. Turnbull, of the Waterways Experiment Station, 
and whoever turns out to be the Chairman of the Soil Mechanics Group of the 
Texas Section, ASCE, at that time. A hearty welcome and an interesting meet- 
ing await all who can attend. 


KNOXVILLE CONVENTION 


University of Tennessee June 4-8, 1956 


E. A. Whitehurst is in charge of the planning of the soil mechanics program. 
He can be reached at 112 Perkins Hall, University of Tennessee, Knoxville, 
Tennessee. 


APRIL NEWSLETTER 


Deadline date for arrival at this office of contributions for the April, 1956, 
Newsletter: February 20, 1956, please. 


Howard P. Hall, Editor 


Department of Civil Engineering = 
Northwestern University 
Evanston, Illinois 
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DIVISION ACTIVITIES 
SOIL MECHANICS AND FOUNDATIONS DIVISION 
Proceedings of the American Society of Civil Engineers 


Soil Mechanics and Foundations Division 
Committee on Publications 


INSTRUCTIONS FOR AUTHORS 
Effective date: 27 August 1955 


The following abbreviated interpretation of requirements has been com- 
piled by the Division Committee on Publications to assist authors in the prep- 
aration of papers for the Soil Mechanics and Foundations Division. It is based 
on Proceedings Separate 290, “Publication Procedure for Technical Papers,” 
as amended by a number of subsequent instructions and changes. 

Papers selected by the Soil Mechanics and Foundations Division for publi- 
cation (In its quarterly journal) constitute a part of the Proceedings of the 
Society. As such, they provide a forum on subjects of current discussion and 
interest. 

While these instructions apply only to preparation of papers for the Pro- 
ceedings, it should be kept in mind that all Transactions papers must originate 
in the Proceedings. Periodically, a very limited number of papers are select- 
ed by the Publications Committee of the Society Board of Direction to provide 
a permanent record of those papers of outstanding and lasting quality. Papers 
potentially in this class should therefore be prepared originally to meet the 
standards of Transactions. These standards are explained in Part II of 
Proceedings-Separate 290. 


General Requirements 


Papers selected for Proceedings should (1) be reasonably clear and con- 
cise, (2) meet minimum requirements of form, and (3) be of interest and con- 
structive value to civil engineers, primarily the members of the Soil Mechan- 
ics and Foundations Division. They need not fulfill Transactions requirements 
that new or original material be presented. However, they should avoid per- 
sonal references, commercial implications, and extensive reprinting of mate- 
rial in print elsewhere. 


Requirements for Mechanical Preparation 


All papers presented should require a minimum of technical editing since 
those selected for Proceedings are returned to headquarters for re-typing of 


Note: No. 1956-3 is part of the copyrighted Journal of the Soil Mechanics and Founda- 
tions Division of the American Society of Civil Engineers, Vol. 82, SM 1, January, 
1956. 
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text, and photo-offset printing, without further editorial correction. Lacking 


an editorial staff, the Division Committee on Publications may elect to return 
without review those papers which are grossly deficient in meeting the pre- 
scribed mechanical requirements. 


A summary of these requirements is presented below. 


Employment. 


Abstract ... 


Figures.... 


Tables. 


Mathematics... 


usually not exceeding 12,000 word-equivalents (22 pages 
as might be printed in Transactions.) Special approval by 
the Society required for papers longer than 18,000 words. 


Not longer than 50 letters and spaces. 


footnote at bottom of first page showing author’s employ- 
ment. 


approximately 50 words, furnished by the author, for use 
in Civil Engineering. 


double-spaced original copy or clear reproduction on 
8-1/2 x 11 in. or 8 x 10-1/2 in. paper (will be retyped by 
printer). 


will be printed with 100 to 69 reduction directly from 
author’s copy; although not desirable, other reduction 
possible if uniform for all figures. 


Black on white background, within an invisible 6-1/2 x 
10-1/2 in. frame. 


Lettering minimized, all descriptive material in text, 
lettered capitals 3/16 in. high, all others in proportion. 


Number of illustrations, especially photographs, not ex- 
cessive; each adequately explained in text. 


Arabic numbering, drawings consecutive with photographs. 


will be printed with 100 to 69 reduction directly from 
author’s copy. 


Black on white background, within invisible frame 6-1/2 x 
10-1/2 in. 


Typing clear and black, suitable for photo reproduction 
Arabic numbering 


will be printed with 100 to 69 reduction directly from 
author’s copy 


Black on white paper, lettered capitals 3/16 in. high, 
other symbols in proportion, or typed as required for 
tables. 


Maximum length of each line 6-1/2 in. 


Procedure for Submitting Individual Papers 

Three copies of the prepared paper should be sent to the Manager of Tech- 
nical Publications, American Society of Civil Engineers, 33 West 39th Street, 
New York 18. They should be accompanied by a letter of transmittal directing 
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1956-3 Soil Mechanics and Foundations Division 
the paper to the Soil Mechanics and Foundations Division and stating that 


a) the manuscript has been presented for the exclusive use of the Society, 
b) it has not been printed, in whole or in part, anywhere else, and 
c) it may be released to the Society’s printer for re-typing and printing. 


If the paper is presented orally, the place and date should be indicated. 


Procedure for Submitting Groups of Papers 


Related papers, all of which are under the jurisdiction of this Division, can 
be published as consecutive papers in the Proceedings if the following pro- 
visions are met. 


a) The authors are in agreement that they should be so published. 

b) A coordinator is selected by those concerned to collect the papers, 
transmit them to the Manager of Technical Publications as a group, 
and specifically request that they be so handled. A copy of the letter of 
transmittal should be sent directly to the Chairman of the Division Pub- 
lication Committee. 

c) The coordinator should supply an over-all explanatory statement of 
about 250 or 300 words to be duplicated on the front page of each pub- 
lished paper. 


Each paper must meet all of the requirements for individual publication. 
Only those papers approved by the Committee on Publications will be included 
in the published group. 


Discussions of Published Papers 

A discussion period of three months from the first date of issue of a paper 
will be allowed. 

Discussers should send the original and one carbon copy of their contribu- 
tions directly to the Manager of Technical Publications. Standards for dis- 
cussion copy will be the same as for papers. 

Each discussion must be directed to a single paper. Where it is desired to 
discuss more than one paper in a group of related Proceedings papers, sepa- 
rate discussions are necessary. 

A discussion should be confined to the scope set by the author in his paper; 
and the discusser should discuss rather than extend the author’s material, 
avoiding side issues and connecting his comments by specific references to 
the main paper. 
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PROCEEDINGS PAPERS 


The technical papers published in the past year are identifted by number below. Tech- 
nical-division sponsorship is indicated by an abbreviation at the end of each Paper Number, 
the symbols referring to: Air Transport (AT), City Planning (CP), Construction (CO), 
Engineering Mechanics (EM), Highway (HW), Hydraulics (HY), Irrigation and Drainage (IR), 
Power (PO), Sanitary Engineering (SA), Soil Mechanics and Foundations (8M), Structural 
(ST), Surveying and Mapping (SU), and Waterways (WW) divisions. Papers sponsored by 
the Board of Direction are identified by the symbols (BD). For titles and order coupons, 
refer to the appropriate issue of “Civil Engineering.” Beginning with Volume 82 (January 
1956) papers were published in Journals of the various Technical Divisions. To locate 
papers in the Journals, the symbols after the paper numbers are followed by a numeral 
designating the issue of a particular Journal in which the paper appeared. For example, 
Paper 861 is identified as 861 (SM1) which indicates that the paper is contained in issue 1 
of the Journal of the Soil Mechanics and Foundations Division. 


VOLUME 81 (1955) 


JANUARY: 583(ST), 584(ST), 585(ST), 586(ST), 587(ST), 588(ST), 589(ST)°, 590(SA), 591(SA), 
592(SA), 593(SA), 594(SA), 595(SA)°, 596(HW), 597(HW), 598(HW)°,599(CP), 600(CP), 601(CP), 
602(CP), 603(CP), 604(EM), 605(EM), 606(EM)°, 607(EM). 


FEBRUARY: 608(WW), 609(WW), 610(WW), 611(WW), 612(WW), 613(WW), 614(WW), 615(Ww), 
€16(WW), 6170R), 618(IR), 619GR), 620(R), 6210R)°, 622(1R), 62300R), 624(HY)*, 625(HY), 
626(HY), 627(HY), 628(HY), 629(HY), 630(HY), 631(HY), 632(CO), 633(CO). 


MARCH: 634(PO), 635(PO), 636(P0), 637(PO), 638(PO), 639(PO), 640(PO), 641(PO)", 642(8A), 
643(SA), 644(SA), 645(SA), G46(SA), G647(SA)°, $48(ST), 649(ST), 650(ST), 651(ST), 652(ST), 
653(ST), 654(ST)°, 655(SA), 656(SM)°, 657(SM;°, 658(SM)°. 


APRIL: 659(ST), 660(ST), 661(ST)©, 662(ST), 663{ST), 664(ST)°, 665(HY)°, 666(HY), 667(HY), 
v68(HY), 669(HY), 670(EM), 671(EM), 672(EM), 673(EM), 674(EM), 675(EM), 676(EM), 677(EM), 
678(HY). 


MAY: 679(ST), 680(ST), 681(ST), 682(ST)°, 683(ST), 684(ST), 685(SA), 686(SA), 687(SA), 688(SA), 
689(SA)°, 690(EM), 691(EM), 692(EM), 693(EM), 694(EM), 695(EM). 696(PO), 697(PO), 698(SA), 
699(PO)¢, 700(PO), 701(ST)*. 


TUNE: 702(HW), 703(HW), 704¢HW)°, 705(IR), 706(IR), 707(IR), 708(IR), 709(HY)*, 710(CP), 
T11(CP), 712(CP), 713(CP)°, 714(HY), 715(HY), 716(HY), 717(HY), 718(SM)°, 719(HY)S, 
720(AT), 721(AT), 722(SU), 723(WW), 724(Ww), 725(ww), 726(ww)°, 727(WW), 728(IR), 
729(IR), 730(SU)°, 731(SU), 


JULY: 732(ST), 733(ST), 734(ST), 735(ST), 736(ST), 737(PO), 738(PO), 739(PO), 740(PO), 
'141(PO), 742(PO), 743(HY), 744(HY), 745(HY), 746(HY), 747(HY), 748(HY)°, 749(SA), 750(SA), 
7S1(SA), 752(8A)°, 753(SM), 754(SM), 755(SM), 756(SM), 757(SM), 758(CO)°, 759(SM)*, 
160(WW)”, 


AUGUST: 1761(BD), 762(ST), 763(ST), 764(ST), 765(ST)©, 766(CP), 767(CP), 76&CP), 769(CP), 
T70(CP), TT1(EM), 772(EM), 3A), 774(EM), T75(EM), 776(EM)°, 777(AT), T78(AT), 
TI9(SA), 780(SA), 781(SA), 782(Sa)<, 783(HW), 784(HW), 785(CP), 786(ST). 


SEPTEMBER: 787(PO), 788(IR), 789(HY), 790(HY), 791(HY), 792(HY), 793(HY), 794(HY)°, 
195(EM), 796(EM), 797(EM), 798{EM), 799(EM)°, 800(WW), 801(WW), 802(WW), 803(WW), 
304(WW), 805(WW), 806(HY), 807(PO)°, 808(IR)°. 


OCTOBER: 809 (ST), 810 (HW)°, 811(ST), 812 (ST)*, 813 (ST)©, 814(EM), 815(EM), 816(EM), 
317(EM), 818(EM), 819(EM)°, 820(SA), 821(SA), 822(SA)°, 823(HW), 824(HW). 


NOVEMBER: 825(3T), 826(HY), 827(ST), 828(ST), 829(sT), 330(sT), 831(sT)°, 832(CP), 
a33(CP), 834(CP), 835(CP)®, 836(HY), 837(HY), 838(HY), 839(HY), 840(HY), 


DECEMBER: 842(SM), 843(SM)°, 844(SU), 845(SU)©, 846(SA), 847(SA), 848(SA)°, 849(ST)°, 
850(ST), 851(ST), 852(ST), 853(ST), 854(CO), 855(CO), 856(CO)°, 857(SU), 858(BD), 859(BD), 
860(BD). 

VOLUME 82 (1956) 


JANUARY: 861(SM1), 862(SM1), 663(EM1), 864(SM1), 865(SM1), 866(5M1), 867(SM1), 868 
(HW1), 869(ST1), 870(EM1), 871(HW1), 872(HW1), 87(HW1), 874(HW1), 875(HW1), 876 
(EM1)°, 877(HWw1)°, 878(sT1)°. 


e. Discussion of several papers, grouped by Divisions. 
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